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14.  ABSTRACT 

The  objective  of  this  thesis  was  to  determine  the  feasibility  of  using  loop  heat  pipes  to  dissipate  waste  heat  from  power 
electronics  to  the  skin  of  a  fighter  aircraft  and  examine  the  performance  characteristics  of  a  titanium-water  loop  heat  pipe 
under  stationary  and  elevated  acceleration  fields.  In  the  past,  it  has  been  found  that  the  boundary  condition  at  the 
condenser  can  be  a  controlling  factor  in  the  overall  performance  of  this  type  of  thermal  management  scheme.  Therefore, 
the  heat  transfer  removed  from  the  aircraft  skin  has  been  determined  by  modeling  the  wing  as  a  flat  plate  at  zero- 
incidence  as  a  function  of  the  following  parameters:  airspeed:  0.8  <  Ma*  <  1.4;  altitude:  0  <  H<  22  km;  wall 
temperature:  105  <  Tw  <  135  °C.  In  addition,  the  effects  of  the  variable  properties  of  air  have  been  taken  into  account. 
Heat  transfer  due  to  thermal  radiation  has  been  neglected  in  this  analysis  due  to  the  low  skin  temperatures  and  high 
airspeeds  up  to  Ma*  =  1.4.  It  was  observed  that  flight  speed  and  altitude  have  a  significant  effect  on  the  heat  transfer 
abilities  from  the  skin  to  ambient,  with  heat  rejection  becoming  more  difficult  with  increasing  Mach  number  or 
decreasing  altitude.  Abstract  concludes  on  reverse 
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14.  ABSTRACT  (concluded) 


An  experiment  has  been  developed  to  examine  operating  characteristics  of  a  titanium-water  loop  heat  pipe  (LHP)  under 
stationary  and  elevated  acceleration  fields.  The  LHP  was  mounted  on  a  2.44  m  diameter  centrifuge  table  on  edge  with  heat 
applied  to  the  evaporator  via  a  mica  heater  and  heat  rejected  using  a  high-temperature  polyalphaolefm  coolant  loop.  The 
LHP  was  tested  under  the  following  parametric  ranges:  heat  load  at  the  evaporator:  100  <  gin  <  600  W;  heat  load  at  the 
compensation  chamber:  0  <  gcc  <  50  W;  radial  acceleration:  0  <  Or  <  10  g.  For  stationary  operation  {a^=  1.0  g,  flr  =  0  g), 
the  LHP  evaporative  heat  transfer  coefficient  decreased  monotonically,  thermal  resistance  decreased  to  a  minimum  then 
increased,  and  wall  superheat  increased  monotonically.  Heat  input  to  the  compensation  chamber  was  found  to  increase  the 
evaporative  heat  transfer  coefficient  and  decrease  thermal  resistance  for  gin  =  500  W.  Flow  reversal  in  the  LHP  was  found 
for  some  cases,  which  was  likely  due  to  vapor  bubble  formation  in  the  primary  wick.  Operating  the  LHP  in  an  elevated 
acceleration  environment  {a^  =  1.0  g,  Or  >  0  g)  revealed  dry-out  conditions  from  gin  =  100  to  400  W  and  varying 
accelerations  and  the  ability  for  the  LHP  to  reprime  after  an  acceleration  event  that  induced  dry-out.  Evaporative  heat 
transfer  coefficient  and  thermal  resistance  as  found  not  to  be  significantly  dependent  on  radial  acceleration.  However,  wall 
superheat  was  found  to  increase  slightly  with  radial  acceleration. 
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1.  CONVECTIVE  HEAT  TRANSFER  FROM  HIGH-SPEED  AIRCRAFT  SKIN 


1.1.  Abstract 

The  objective  of  the  present  analysis  was  to  determine  the  feasibility  of  using  loop 
heat  pipes  to  dissipate  waste  heat  from  power  electronics  to  the  skin  of  a  fighter  aircraft. 
In  the  past,  it  has  been  found  that  the  boundary  condition  at  the  condenser  can  be  a 
controlling  factor  in  the  overall  performance  of  this  type  of  thermal  management  scheme. 
Therefore,  the  heat  transfer  removed  from  the  aircraft  skin  has  been  determined  by 
modeling  the  wing  as  a  flat  plate  at  zero-incidence  as  a  function  of  the  following 
parameters:  airspeed:  0.8  <  Maoo  <  1.4;  altitude:  0  <//<  22  km;  wall  temperature:  105 
<  Ik  <  135°C.  In  addition,  the  effects  of  the  variable  properties  of  air  have  been  taken 
into  account.  Heat  transfer  due  to  thermal  radiation  has  been  neglected  in  this  analysis 
due  to  the  low  skin  temperatures  and  high  airspeeds  up  to  Maoo  =  1.4.  It  was  observed 
that  flight  speed  and  altitude  have  a  significant  effect  on  the  heat  transfer  abilities  from 
the  skin  to  ambient,  with  heat  rejection  becoming  more  difficult  with  increasing  Mach 
number  or  decreasing  altitude. 

1.2.  Introduction 

The  More  Electric  Aircraft  initiative  (MEA)  is  the  concept  for  future  aircraft 
including  warfighter,  transport,  helicopters,  and  commercial  aircraft.  This  approach  has 
been  adopted  by  the  United  States  Air  Eorce  since  the  early  1990’s  with  the  purpose  of 
reducing  or  removing  as  many  of  the  hydraulic,  mechanical,  and  pneumatic  power 
components  and  replacing  them  with  electrically  driven  devices.  This  approach  to 
aircraft  design  was  first  envisioned  during  World  War  II.  However,  at  that  time,  the 
power  generation  capability  and  power  conditioning  equipment  required  was  not  feasible 
due  to  volume  requirements.  As  a  result,  hydraulic,  pneumatic,  and  mechanical  systems 
became  the  norm  for  aircraft  until  this  initiative.  Under  the  MEA  paradigm,  power  for 
systems  such  as  flight  control  actuation,  anti-ice,  braking,  environmental  control,  engine 
starting,  and  fuel  pumping  will  be  provided  by  a  starter/generator  driven  by  the  gas 
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generator  spool  of  the  aireraft  engine  (Quigley,  1993).  The  MEA  initiative  has  been 
analytieally  proven  to  improve  aireraft  reliability,  maintainability,  support,  and  operations 
eost  as  well  as  reduee  weight,  volume,  and  enhanee  battle  damage  reeonfigurability 
(Cloyd,  1997). 

While  the  reduetion  of  hydraulie,  pneumatie,  and  meehanieal  systems  in  favor  of 
eleotrieal  systems  is  benefieial,  it  presents  a  problem  in  terms  of  thermal  management. 
Replaeing  the  eentralized  hydraulie  system  with  an  eleetrieal  based  system  removes  a 
primary  method  of  transporting  and  removing  waste  heat  (Vrable  and  Yerkes,  1998).  A 
separate  eooling  fluid  system  for  thermal  management  would  be  eontrary  to  the  goals  of 
the  MEA  initiative.  Therefore,  thermal  management  would  need  to  be  distributed  over 
the  entire  aireraft.  As  a  result,  a  new  approaeh  to  thermal  management  involves  handling 
heat  loads  on  a  loeal  level.  This  means  taking  individual  eomponents  in  the  aireraft  and 
loeally  handling  their  heat  rejeetion  requirements. 

The  operating  envelope  for  military  aireraft  plaees  stringent  limitations  on  any 
proposed  thermal  management  system.  The  on-board  eleetrieal  flight  eontrol  aetuation 
system  operates  at  altitudes  from  sea  level  to  above  12  km,  airspeeds  from  stationary  to 
supersonic  speeds,  transient  body  forces  up  to  9  g  due  to  maneuvering,  and  ambient 
temperatures  from  -68  to  58°C.  MEA  has  resulted  in  the  development  of  high- 
temperature,  high-efficiency,  and  high-density  power  electronic  component  technologies. 
The  next-generation  power  electronics  will  be  capable  of  operating  at  cold  plate 
temperature  excursions  up  to  200“C,  which  presents  an  opportunity  to  reject  heat  through 
the  aircraft  skin  to  the  ambient  using  passive  cooling.  In  addition,  the  actuation  system 
rejects  heat  continuously  at  a  rate  of  g  =  500  W  (^  =  3  W/cm  )  and  has  transient  heat 
rejection  rates  of  Q  =  5000  W  over  a  period  of  one  second.  Possible  thermal 
management  scenarios  include  direct  connection  of  the  electronics  package  to  the  skin, 
high-thermal  conductivity  graphite  straps,  or  the  use  of  a  loop  heat  pipe  between  the 
package  and  the  skin  to  provide  mounting  flexibility.  The  objective  of  this  analysis  is  to 
determine  the  external  heat  transfer  possibilities  of  the  aircraft  skin.  The  heat  flux  and 
heat  transfer  coefficient  have  been  found  as  functions  of  the  skin  and  ambient 
temperatures,  the  altitude,  and  airspeed. 
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1,3,  Mathematical  Model 

The  temperature  and  density  of  air  vary  eonsiderably  with  altitude  and  also  vary 
day-to-day  depending  on  weather  eonditions.  In  order  to  be  eonservative  in  the 
ealculation  of  heat  transfer  coefficients,  data  for  the  highest  temperature  recorded  with  a 
frequency-of-occurrence  of  1%  were  used  to  generate  equations  for  temperature  and 
density  versus  altitude  (DOD,  1997)  as  shown  in  Figure  1.1  and  Table  1.1.  Also 
presented  are  data  for  the  lowest  temperature  recorded  with  a  frequency-of-occurrence  of 
1%  (DOD,  1997)  and  data  for  the  “standard  atmosphere”  (Anderson,  2000). 

The  film  temperature  was  used  as  the  reference  temperature  to  evaluate  the  air 
properties  (White,  1988) 

r  =  r,o(0.5  +  0.039Mai)  +  0.5r^  (1.1) 

The  air  density  at  the  film  temperature  and  at  altitude  was  evaluated  using  the  perfect  gas 
law 


P‘=P.»(p)  (1-2) 

The  freestream  speed  of  sound  is 

ttoo  =  ^yRTa,  (1-3) 

The  freestream  velocity  is 


fioo  =  Maoottoo 


(1.4) 


The  absolute  viscosity  of  air  is  given  by  the  following  relation  (NACA,  1953) 

0.76 


M  =  Mr 


4) 


(1.5) 


where  //r  is  a  reference  viscosity  evaluated  at  a  known  reference  temperature  7r. 

The  Reynolds  number  for  a  plate  of  length  L  is  determined  by  evaluating  the 
properties  of  air  at  the  freestream  condition. 

Poo  UqqL 


Rci  — 


(1.6) 


3 


Regression  equations  for  the  speeifie  heat  and  Prandtl  number  were  determined  as 
funetions  of  temperature  using  data  from  Ineropera  and  DeWitt  (2002),  as  shown  in  Table 

1.2. 

The  adiabatie  wall  temperature  is  (White,  1988) 


T  =T 

aw  oc 


1  +  r 


(^) 


Ma? 


where  the  reeovery  faetor  is 


(1.7) 


r  —  {  for  laminar  flow  ^  g^ 

[pri/3  fQj,  turbulent  flow 

For  the  purposes  of  this  analysis,  Reynolds  numbers  less  than  500,000  were  eonsidered  to 
be  laminar,  greater  than  500,000  were  turbulent.  The  loeal  skin  friction  coefficient  at  the 
end  of  the  plate  was  found  by  evaluating  the  air  properties  at  the  film  temperature.  For 
laminar  flow,  the  skin  friction  coefficient  is  given  by  (White,  1988) 


0.664 


(1.9) 


and  for  turbulent  flow 


0.455 

^0.06p*l/^Lj  (1.10) 

The  local  Stanton  number  at  the  end  of  the  plate  for  laminar  flow  is  given  by  (White, 
1988) 

StL  =  0.332ReL^^Pr-2/^  (1.11) 

and  for  turbulent  flow 


K  _ c|y2 _ 

1  +  12.7(Pr’2^3  _ 

The  local  heat  transfer  coefficient  at  the  end  of  the  plate  is 


/iL  =  Stip*u^c; 


(1.12) 


(1.13) 
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The  local  heat  transfer  coefficient  was  calculated  using  the  appropriate  skin  friction 
coefficient  and  Stanton  number  based  on  laminar  or  turbulent  flow.  The  average  heat 
transfer  coefficient  over  the  length  of  the  plate  is  approximated  by  (White,  1988) 

h  =  1.15/iL  (1.14) 

The  heat  flux  dissipated  over  the  plate,  both  local  and  average,  is  defined  in  terms  of  the 
adiabatic  wall  temperature  (White,  1988) 

qw  =  KTw-T^^)  (1-15) 

Thermal  radiation  was  neglected  in  this  analysis  as  it  contributed  less  than  1.6%  to  the 
total  heat  rejected  from  the  plate  surface. 


1.4.  Results  and  Discussion 

The  adiabatic  wall  temperature  is  shown  in  Figure  1.2  as  a  function  of  altitude  and 
Mach  number.  The  overall  trend  of  the  adiabatic  wall  temperature  with  altitude  follows 
the  freestream  air  temperature  in  Figure  1.1  and  increases  with  Mach  number  as 
expected.  Figure  1.3  presents  the  temperature  difference,  AT=  {Taw  -  Too),  versus  altitude. 
This  temperature  difference  demonstrates  the  increase  in  the  adiabatic  wall  temperature 
over  the  freestream  due  to  aerodynamic  heating.  The  temperature  difference  AT  =  ( Jk  - 
Taw)  is  given  in  Figure  1.4.  Of  interest  is  the  portion  of  the  curves  in  which  this 
difference  is  negative,  which  indicates  that  heat  is  transferred  from  the  air  to  the  aircraft 
skin.  The  maximum  Mach  number  achievable  before  heat  is  transferred  from  the  air  to 
the  skin  is  given  by 


Ma 


00, max 


(1.16) 


and  is  plotted  in  Figure  1.5  over  a  range  of  wall  temperatures.  The  maximum  Mach 
number  increases  with  altitude  and  wall  temperature  up  to  a  maximum  at  approximately 
18  km.  In  Figure  1.6,  the  average  convective  heat  transfer  coefficient  decreases 
monotonically  with  altitude  due  to  the  continual  decrease  in  the  air  density.  In  general, 
the  convective  heat  transfer  coefficient  increases  with  Mach  number,  as  expected.  The 
average  heat  flux  dissipated  from  the  plate  is  shown  in  Figure  1.7.  For  low  Mach 
numbers,  the  heat  flux  is  positive  for  all  values  of  altitude,  which  indicates  that  heat  is 
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transferred  from  the  aircraft  skin  to  the  air.  At  high  Mach  numbers,  however,  the  heat 
flux  is  negative  at  low  altitudes  due  to  the  negative  Ar  as  shown  in  Figure  1.4.  This 
means  that  the  adiabatic  wall  temperature  is  higher  than  the  skin  temperature  due  to 
aerodynamic  heating  effects.  The  effect  of  heated  plate  length  on  the  local  heat  flux  for 
//=  0,  10,  and  20  km  is  shown  in  Figure  1.8.  The  local  heat  flux  starts  low  and  decreases 
in  the  laminar  region  of  the  plate,  and  then  increases  as  the  flow  transitions  to  turbulent 
where  it  once  again  decreases.  In  general,  the  average  heat  flux  follows  the  behavior  of 
the  local  heat  transfer  coefficient,  where  Al  is  high  at  the  leading  edge  and  at  the 
beginning  of  turbulent  flow  and  decreases  as  the  boundary  layer  grows.  One  item  to  note 
is  that  Figure  1.9  shows  the  average  heat  flux  dissipated  over  the  plate  versus  altitude  for 
the  1%  hot  day,  the  1%  cold  day,  and  the  standard  atmosphere  data  as  presented  in  Figure 
1.1.  At  low  altitudes,  is  significantly  higher  for  the  1%  cold  day  due  to  the  combined 

effects  of  the  lower  atmospheric  temperature  and  the  higher  air  density.  The  effect  of 
wall  temperature  on  average  heat  flux  for  a  given  airspeed  is  shown  in  Figure  1.10.  The 
heat  flux  increases  dramatically  with  altitude  and  wall  temperature  for  low  altitudes. 

1,5,  Conclusions 

An  analysis  of  the  heat  transfer  from  a  heated  plate  has  provided  important 
insights  for  the  possible  use  of  the  aircraft  skin  to  reject  heat  from  electric  actuator 
systems.  It  was  found  that  the  altitude  and  speed  of  the  aircraft  significantly  affected  the 
amount  of  heat  that  could  be  rejected  from  the  skin.  Aerodynamic  heating  of  the  skin 
reduced  the  heat  transfer,  and  if  the  Mach  number  was  high  enough,  heat  transfer  from 
the  skin  to  the  air  went  to  zero.  A  performance  map  of  this  phenomenon  was  provided. 
The  altitude  of  the  aircraft  affected  the  freestream  temperature  and  density,  which  in  turn 
affected  the  overall  heat  transfer  coefficient.  It  was  also  shown  that  the  assumption  of  a 
“standard  atmosphere”  could  result  in  significant  errors  in  the  prediction  of  the  heat 
dissipation  as  compared  to  the  data  for  the  1%  hot  day  or  the  1%  cold  day.  The  analysis 
showed  that  the  aircraft  skin  temperature,  which  is  directly  influenced  by  the  actuator 
thermal  management  system,  has  a  strong  effect  on  the  heat  dissipation  rate,  especially  at 
low  altitudes. 
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0  5  10  15  20 

H  (km) 

Figure  1.1.  Comparison  of  atmospheric  properties  versus  altitude:  (a)  Temperature;  (b) 
Density  (DOD,  1997;  Anderson,  2000). 
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Figure  1.2.  Adiabatic  wall  temperature  versus  altitude  for  various  Mach  numbers  (1% 
hot  day). 


Figure  1.3.  Temperature  difference  -  T^)  versus  altitude  for  various  Mach  numbers 
(1%  hot  day). 
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Figure  1.4.  Temperature  difference  (T^  versus  altitude  for  various  Mach  numbers 
(rw=  135°C,  1%  hot  day). 
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Figure  1.5.  Maximum  Mach  number  before  heat  is  transferred  from  the  air  to  the  skin 
versus  altitude  for  various  wall  temperatures  (1%  hot  day). 
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Figure  1.6.  Average  convective  heat  transfer  coefficient  versus  altitude  for  various  Mach 
numbers  (JW  =  135“C,  L=  1.0  m,  1%  hot  day). 


Figure  1.7.  Average  heat  flux  dissipated  over  the  plate  versus  altitude  for  various  Mach 
numbers  =  135“C,  L=  1.0  m,  1%  hot  day). 
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Figure  1.9.  Average  heat  flux  dissipated  over  the  plate  versus  altitude  for  various 
atmospheric  conditions  (JW  =  135°C,  L  =  1.0  m,  Maoo  =  0.98)  (DOD,  1997;  Anderson, 
2000). 


Figure  1.10.  Average  heat  flux  dissipation  versus  altitude  for  various  wall  temperatures 
(Maoo  =  0.98,  Z  =  1.0  m,  1%  hot  day). 
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Table  1.1.  Regression  equations  for  air  properties  versus  altitude  for  1%  hot  (DOD, 
1997). 


y  =  ao  +  a\H  + 

{H  in  km) 


Property 

ao 

a\ 

a2 

<33 

a4 

Too  rc) 

4.8507E+1 

-9.5033E+0 

5.3483E-1 

-2.8994E-2 

7.7664E-4 

0.99779 

yOoo  (kg/m^) 

1.0868E+0 

-8.9917E-2 

2.0898E-3 

-4.9336E-6 

— 

0.99954 

Table  1.2.  Regression  equations  for  air  properties  versus  temperature  (Ineropera  and 
DeWitt,  2002). 


_y  =  ao  +  a\T  +  a27^  +  a^'f' 
(TinK) 


Property 

ao 

a\ 

a2 

A3 

Cp  (J/kg-K) 

1.0187E+3 

-6.9921E-2 

-3.3333E-5 

4.4444E-7 

0.99916 

Pr 

8.6418E-1 

-9.4177E-4 

1.7778E-6 

-1.2593E-9 

0.99725 
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2.  TITANIUM-WATER  LOOP  HEAT  PIPE  CHARACTERISTICS  UNDER 
STATIONARY  AND  ELEVATED  ACCELERATION  FIELDS 

2.1,  Abstract 

An  experiment  has  been  developed  to  examine  operating  eharaeteristies  of  a 
titanium-water  loop  heat  pipe  (LHP)  under  stationary  and  elevated  acceleration  fields. 
The  LHP  was  mounted  on  a  2.44  m  diameter  centrifuge  table  on  edge  with  heat  applied 
to  the  evaporator  via  a  mica  heater  and  heat  rejected  using  a  high- temperature 
polyalphaolefin  coolant  loop.  The  LHP  was  tested  under  the  following  parametric 
ranges:  heat  load  at  the  evaporator:  100  <  Q\n  <  600  W;  heat  load  at  the  compensation 
chamber:  0  <  gcc  <  50  W;  radial  acceleration:  0  <  ar  <  10  g.  For  stationary  operation 
=  1.0  g,  ar  =  0  g),  the  LHP  evaporative  heat  transfer  coefficient  decreased  monotonically, 
thermal  resistance  decreased  to  a  minimum  then  increased,  and  wall  superheat  increased 
monotonically.  Heat  input  to  the  compensation  chamber  was  found  to  increase  the 
evaporative  heat  transfer  coefficient  and  decrease  thermal  resistance  for  =  500  W. 
Flow  reversal  in  the  LHP  was  found  for  some  cases,  which  was  likely  due  to  vapor 
bubble  formation  in  the  primary  wick.  Operating  the  LHP  in  an  elevated  acceleration 
environment  revealed  dry-out  conditions  from  =  100  to  400  W  and  varying 
accelerations  and  the  ability  for  the  LHP  to  reprime  after  an  acceleration  event  that 
induced  dry-out.  Evaporative  heat  transfer  coefficient  and  thermal  resistance  was  found 
not  to  be  significantly  dependent  on  radial  acceleration.  However,  wall  superheat  was 
found  to  increase  slightly  with  radial  acceleration. 

2.2.  Introduction 

Loop  heat  pipes  (LHP's)  are  two-phase  thermal  transport  devices  that  operate 
passively  using  the  latent  heat  of  vaporization  to  transport  heat  from  one  location  to 
another.  The  LHP  was  invented  in  1972  by  Gerasimov  and  Maidanik  (Maidanik,  2005) 
in  the  former  Soviet  Union,  and  was  later  patented  in  the  United  States  (Maidanik  et  ah, 
1985).  The  LHP  consists  of  an  evaporator,  compensation  chamber,  liquid  and  vapor 
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transport  lines  made  of  smooth  tubing,  and  a  condenser  as  shown  in  Figure  2.1.  Heat  is 
applied  directly  to  the  exterior  wall  of  the  evaporator,  which  often  has  a  circular  cross- 
section.  The  majority  of  the  input  heat  is  used  to  vaporize  the  working  fluid  within  the 
primary  wick  structure,  which  is  an  inverted  meniscus  wick  in  direct  contact  with  the 
exterior  evaporator  wall.  The  vapor  is  captured  in  the  axial  vapor  grooves  in  the  primary 
wick  and  is  directed  via  a  manifold  at  the  end  of  the  evaporator  to  the  vapor  line  due  to 
the  increased  pressure  within  the  evaporator.  Due  to  evaporation,  menisci  are  developed 
in  the  primary  wick  which  establishes  a  capillary  pressure  head  that  returns  liquid  to  the 
evaporator  from  the  condenser.  This  capillary  head  must  be  greater  than  the  total  system 
pressure  drop  in  order  for  the  LHP  to  continue  to  operate  without  drying  out. 

The  vapor  from  the  evaporator  section  travels  via  the  vapor  line  to  the  condenser 
section,  which  is  also  made  of  smooth  tubing.  Heat  is  rejected  from  the  condenser  to  the 
ultimate  heat  sink.  The  working  fluid  enters  the  condenser  as  a  superheated  vapor.  After 
sufficient  heat  is  rejected,  the  vapor  becomes  a  saturated  vapor,  a  two-phase  mixture,  a 
saturated  liquid,  and,  depending  on  the  amount  of  heat  rejection,  it  may  or  may  not 
become  a  subcooled  liquid.  The  location  of  the  point  at  which  the  working  fluid  becomes 
a  subcooled  liquid  (2(p-l(p)  is  dependent  on  the  heat  input  at  the  evaporator,  the  heat 
rejection  at  the  condenser,  and  the  saturation  temperature  in  the  compensation  chamber. 
After  exiting  the  condenser  section,  the  liquid  will  continue  to  lose  heat  due  to  convection 
and/or  thermal  radiation  to  the  ambient.  The  subcooled  liquid  returns  to  the  evaporator 
via  the  bayonet  tube,  which  delivers  the  liquid  to  the  end  of  the  evaporator  where  the 
vapor  manifold  resides. 

As  stated  previously,  most  of  the  evaporator  heat  input  evaporates  liquid  in  the 
primary  wick.  The  rest  of  the  heat  is  transferred  by  conduction  through  the  primary  wick, 
where  liquid  is  evaporated  into  vapor  channels  leading  to  the  compensation  chamber 
(Figure  2.2).  Part  of  this  vapor  stream  condenses  onto  the  secondary  wick,  which  is  in 
intimate  contact  with  the  bayonet  tube.  This  heat  transfer  to  the  bayonet  tube  raises  the 
temperature  of  the  subcooled  liquid  entering  the  compensation  chamber  to  the  saturation 
temperature  as  it  travels  to  the  end  of  the  evaporator.  The  rest  of  the  vapor  condenses 
onto  the  wick  lining  the  compensation  chamber.  This  latent  heat  is  then  rejected  from  the 
compensation  chamber  to  the  ambient.  The  condensate  in  the  compensation  chamber  is 
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drawn  back  to  the  evaporator  seetion  through  the  seeondary  wiek  by  eapillary  aetion.  In 
this  way,  the  seeondary  wiek  and  the  eompensation  ehamber  behave  similar  to  a 
eonventional  heat  pipe. 

The  eompensation  ehamber  allows  the  LHP  to  automatieally  regulate  itself  during 
transient  situations  like  startup,  shutdown,  or  a  ehange  in  the  operating  eonditions.  The 
eompensation  ehamber  provides  for  storage  of  exeess  liquid  when  the  evaporator  heat 
input  is  high,  where  the  majority  of  the  eondenser  seetion  is  free  of  subeooled  liquid.  The 
eompensation  ehamber  ean  also  be  used  to  eontrol  the  loeation  of  the  2(p-l(p  point  in  the 
eondenser.  Controlling  the  heat  transfer  through  the  shell  of  the  eompensation  ehamber 
ean  adjust  the  saturation  point  in  the  eondenser,  thereby  ehanging  the  amount  of 
subeooling  of  the  liquid  returning  to  the  evaporator. 

There  has  been  limited  experimentation  on  the  aeeeleration  effeets  on  loop  heat 
pipes  and  heat  pipes.  Ku  et  al.  (2000a)  performed  experiments  on  a  miniature 
aluminum/anhydrous  ammonia  LHP  by  using  a  spin  table  to  examine  the  effeets  of 
varying  aeeeleration  on  start-up.  Four  mounting  eonfigurations  were  examined:  (1) 
horizontally  with  the  eompensation  ehamber  and  liquid  line  outboard  on  the  table,  (2) 
horizontally  with  the  evaporator  and  vapor  line  outboard  on  the  table,  (3)  vertieally  with 
evaporator  above  the  eompensation  ehamber  with  no  radial  aeeeleration,  and  (4) 
vertieally  with  evaporator  below  the  eompensation  ehamber  with  no  radial  aeeeleration. 
Several  different  experiments  were  eondueted,  ineluding  LHP  startup  before  aeeeleration 
was  applied  and  viee  versa,  as  well  as  varying  heat  load  inputs  up  to  Qm  =  100  W. 
Several  aeeeleration  profdes  were  examined,  ineluding  =  0.0  g,  eonstant  =  1.2  g, 
eonstant  Ur  =  4.8  g,  eombination  of  eonstant  ar  =  1.2  and  4.8  g,  eonstant  ar=  1.2  g  for  30 
seeonds  followed  by  =  0.0  g  for  300  seeonds  periodieally,  eonstant  =  4.8  g  for  30 
seeonds  followed  by  =  0.0  g  for  300  seeonds  periodieally,  and  eombinations  of  =1.2 
and  4.8  g  followed  by  ar  =  0.0  g  for  300  seeonds  periodieally.  Their  experimental  results 
indieated  that  the  wall  superheat,  defined  as  the  differenee  between  the  evaporator  and 
eompensation  ehamber  wall  temperatures,  appeared  to  be  independent  of  input  heat  load 
and  aeeeleration.  When  temperature  overshoot  in  the  evaporator  was  examined,  for  heat 
loads  greater  than  =  50  W,  there  was  essentially  no  overshoot.  For  smaller  heat  loads, 
sueh  as  at  =  5  W,  a  temperature  overshoot  of  a  few  degrees  was  always  observed,  but 
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at  Qin  =  25  W,  the  temperature  overshoot  ranged  from  0  to  45°C.  In  every  experiment, 
the  LHP  started  sueeessfully. 

Ku  et  al.  (2000b),  in  an  extension  of  the  previous  experimental  study,  examined 
the  temperature  stability  of  the  same  miniature  LHP  under  varying  heat  loads  and 
aeeeleration  levels.  Their  experimental  results  showed  that  the  radial  aeeeleration  eaused 
a  redistribution  of  fluid  in  the  evaporator,  eondenser,  and  eompensation  ehamber.  This  in 
turn  ehanged  the  LHP  operating  temperature.  The  effeet  was  not  universal,  in  the  sense 
that  all  the  operating  eonditions  needed  to  be  taken  into  aeeount.  With  suffieient  time, 
eonstant  aeeeleration  eould  either  inerease  or  deerease  the  LHP  operating  temperature. 
Periodie  aeeeleration  led  to  a  quasi-steady  operating  temperature.  Temperature  hysteresis 
eould  also  be  eaused  by  the  radial  aeeeleration.  In  all  of  the  experiments  the  LHP 
eontinued  to  operate  without  problems. 

Similar  researeh  has  been  eondueted  to  examine  body  foree  effeets  on  heat  pipes. 
Ponnappan  et  al.  (1992)  examined  a  flexible  eopper-water  arterial  wiek  heat  pipe 
subjeeted  to  transverse  aeeeleration  using  a  eentrifuge  table.  Evaporator  heat  loads  up  to 
gin  =  150  W  and  steady  state  radial  aeeelerations  up  to  =  10.0  g  were  investigated. 
Transport  eapaeity  of  the  heat  pipe  dropped  from  gout  =  138  W  at  radial  aeeelerations  of 
ar  =  1.0  g  to  gout  =  60  W  at  ar  =  10.0  g.  The  temperature  differenee  between  the 
evaporator  and  eondenser  remained  fairly  eonstant  up  to  ar  =  4  g  then  deereased  from  ar  = 
4  to  10  g.  This  deerease  was  due  to  a  more  uniform  distribution  of  fluid  within  the  wiek 
at  the  higher  radial  aeeeleration. 

Yerkes  and  Beam  (1992)  examined  the  same  flexible  eopper-water  arterial  wiek 
heat  pipe  as  Ponnappan  et  al.  under  transient  transverse  and  axial  aeeeleration  forees  with 
periodie  and  burst  transverse  aeeelerations  from /=  0.01  to  0.03  Hz  and  magnitudes  from 
Or  =  1.1  to  9.8  g  peak-to-peak  and  evaporator  heat  inputs  up  to  gin  =  83  W.  It  was 
observed  that  pooling  of  exeess  fluid  had  a  signifieant  effeet  on  the  heat  transport  of  the 
heat  pipe  at  steady  state  transverse  aeeeleration.  Heat  transport  potential  deereased  with 
inereasing  transverse  aeeeleration  eausing  partial  dry-out  of  the  artery  and  pooling  in  the 
eondenser.  The  heat  pipe  was  able  to  reprime  after  dry-out  events  with  subsequent 
reduetion  of  transverse  aeeeleration.  Under  eyelie  transverse  aeeeleration,  signifieant 
fluid  slosh  was  thought  to  ereate  a  eyelie  variation  in  heat  pipe  temperature.  Temperature 
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rise  was  lower  at  the  onset  of  dry-out  eonditions  when  eompared  to  steady  state 
transverse  aeeeleration.  Frequeney  of  the  steady  periodie  burst  transverse  aeeeleration 
had  no  effeet  on  the  heat  pipe  temperature  and  tended  to  delay  the  onset  of  dry-out. 

Thomas  and  Yerkes  (1996)  examined  the  same  flexible  eopper- water  arterial  wiek 
heat  pipe  as  Ponnappan  et  al.  with  evaporator  heat  loads  from  Q\n  =  75  to  150  W, 
eondenser  temperatures  of  =  3,  20,  and  35°C,  and  sinusoidal  aeeeleration  frequeneies 
of/=  0,  0.01,  0.05,  0.1,  0.15,  and  0.2  Hz.  The  amplitude  of  the  radial  aeeeleration  ranged 
from  ar  =  1.1  to  9.8  g.  The  effeets  of  the  previous  dry-out  history  of  the  heat  pipe  were 
also  examined.  It  was  diseovered  that  the  thermal  resistanee  inereased  and  then 
deereased  with  respeet  to  inereasing  aeeeleration  frequeney.  The  thermal  resistanee  also 
inereased  with  inereasing  evaporator  heat  loads.  The  previous  dry-out  history  adversely 
affeeted  the  thermal  resistanee  of  the  heat  pipe  when  dry-out  oeeurred  prior  to  inereasing 
the  aeeeleration  frequeney. 

Thomas  et  al.  (1998)  examined  a  helieally  grooved  eopper-ethanol  heat  pipe  as  a 
funetion  of  evaporator  heat  input  and  transverse  radial  aeeeleration.  Heat  loads  ranging 
from  gin  =  20  to  250  W  were  applied  to  the  evaporator.  At  gin  =  20  W  the  heat  pipe  did 
not  experienee  any  dry-out  eonditions  when  the  radial  aeeeleration  was  inereased  and 
then  deereased  stepwise  from  ar  =  0  to  10  g.  At  gin  =  50  W,  the  heat  pipe  experieneed 
dry-out  eonditions  at  ar  =  0  and  2  g,  but  quiekly  reprimed  at  the  higher  radial 
aeeelerations.  This  indieated  the  elevated  body  forees  aetually  aided  the  performanee  of 
the  heat  pipe  by  inereasing  the  eapillary  limit  due  to  the  forees  generated  from 
aeeeleration  gradients  down  the  length  of  the  helieal  groove.  The  thermal  resistanee  of 
the  heat  pipe  was  noted  to  deerease  then  inerease  with  inereasing  heat  transported  when 
dry-out  started. 

Zaghdoudi  and  Sarno  (2001)  examined  the  body  foree  effeets  on  a  flat  eopper- 
water  heat  pipe  via  a  eentrifuge  setup.  The  heat  pipe  was  mounted  sueh  that  the 
aeeelerating  forees  were  opposite  to  the  liquid  flow,  or  in  an  “unfavorable”  mounting 
eondition.  Three  types  of  aeeelerations  were  performed  in  this  study:  A  parabolie  profile 
from  ar  =  0  to  10  to  0  g  with  a  5  seeond  stabilization  at  ar  =  10  g,  a  step  inerease  from  = 
0  to  10  to  0  g  with  a  10  seeond  stabilization  at  eaeh  step,  and  inereasing  then  deereasing 
the  aeeeleration  from  ar  =  0  to  10  g  after  thermal  stabilization.  Heat  loads  of  gin  =  20,  40, 
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and  60  W  were  applied  to  examine  the  effect  on  evaporator  and  condenser  temperature  as 
well  as  thermal  resistance.  For  the  first  two  types  of  acceleration  profile,  it  was  observed 
there  was  a  delayed  increase  in  evaporator  temperature  and  decrease  in  condenser 
temperature.  This  was  likely  due  to  the  pooling  of  fluid  in  the  condenser.  Thermal 
resistance  also  experienced  a  delayed  increase  in  onset  and  remained  elevated  even  in  the 
absence  of  an  accelerating  force.  For  the  third  type  of  acceleration  profile,  there  was  a 
much  more  gradual  increase  in  evaporator  temperature  and  nearly  negligible  decrease  in 
condenser  temperature,  quickly  returning  to  normal  in  the  absence  of  the  accelerating 
force.  Thermal  resistance  had  a  similar  trend,  quickly  returning  to  normal  after  the 
acceleration  burst.  This  suggested  that  the  heat  pipe  quickly  reprimed  after  the 
acceleration  event.  These  tests  demonstrated  the  importance  of  prior  operation  history 
when  the  heat  pipe  was  subjected  to  elevated  body  forces. 

The  objective  of  the  present  experiment  was  to  determine  the  operating 
characteristics  of  a  titanium-water  loop  heat  pipe  subjected  to  varying  heat  loads  and 
accelerations.  Transient  temperature  distributions,  the  evaporative  heat  transfer 
coefficient,  and  the  thermal  resistance  have  been  found  in  terms  of  the  heat  input  at  the 
evaporator,  heat  input  at  the  compensation  chamber,  and  radial  acceleration  field.  In 
addition,  the  transient  behavior  during  startup  and  steady  operation  has  been  examined. 
A  performance  map  has  been  developed  that  relates  dry-out  to  the  heat  load  and  radial 
acceleration  for  the  experimental  conditions  described.  The  experimental  parametric 
ranges  were  as  follows:  heat  load  at  the  evaporator:  100  <  Qm  <  600  W;  heat  load  at  the 
compensation  chamber:  0  <  <  50  W;  radial  acceleration:  0  <  ar  <  10  g. 

2,3.  Experimental  Setup 

The  Centrifuge  Table  Test  Bed  at  Wright-Patterson  Air  Force  Base  (AFRL/RZPS) 
was  used  to  determine  the  heat  transfer  characteristics  of  the  titanium-water  LHP  under 
stationary  and  elevated  acceleration  fields.  A  schematic  of  this  test  bed  can  be  seen  in 
Figure  2.3.  The  test  bed  consisted  of  a  2.44  m  diameter  horizontal  rotating  table  driven 
by  a  20  hp  DC  electric  motor.  The  test  bed  was  able  to  deliver  the  following  to  devices 
mounted  to  the  rotating  table:  Conditioned  DC  electrical  power  through  three  separate 
power  supplies,  120  VAC  power,  temperature-controlled  ethylene  glycol  coolant,  and 
electrical  signals  for  analog  or  digital  control.  In  addition,  electrical  signals  were 
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collected  from  instruments  on  the  table  and  stored  in  a  data  acquisition  computer.  The 
radial  acceleration  could  exceed  ar  =  12  g,  with  a  maximum  onset  of  approximately  d^  = 
10  g/s,  inducing  a  tangential  acceleration.  The  acceleration  field  could  be  varied 
manually  using  a  potentiometer,  or  controlled  digitally  using  a  signal  generator  in  the 
data  acquisition  system.  The  acceleration  field  was  measured  using  an  orthogonal  triaxial 
accelerometer  (Columbia  SA-307HPTX)  with  an  uncertainty  of  ±  0.01  g. 

Power  was  supplied  to  heaters  on  the  table  by  three  precision  power  supplies 
(Kepco  ATE150-7M,  Kepco  ATE  1 50-3. 5M,  and  HP  6290A)  through  power  slip  rings. 
These  slip  rings  were  separated  from  the  instrumentation  slip  rings  to  reduce  electrical 
noise.  The  heater  power  was  calculated  by  multiplying  the  voltage  drop  across  the  heater 
by  the  current.  The  current  was  determined  from  the  voltage  drop  across  a  precision 
resistor  in  series  with  the  heater.  This  type  of  measurement  was  required  due  to  the 
voltage  drop  between  the  control  room  and  the  table.  The  uncertainty  in  this 
measurement  was  less  than  2.0% 

Heat  was  rejected  from  the  centrifuge  table  using  an  ethylene -glycol/water 
mixture  that  was  delivered  to  the  rotating  centrifuge  table  via  a  double-pass  hydraulic 
rotary  coupling  (Deublin  1690-000-115).  The  temperature  of  the  coolant  was  maintained 
at  a  constant  setting  by  a  recirculating  chiller  (Neslab  HX-300).  The  volumetric  flow  rate 
of  the  coolant  mixture  was  controlled  using  a  high-pressure  booster  pump,  which  aided 
the  low-pressure  pump  in  the  recirculating  chiller.  Throughout  experimentation  the  flow 
rate  was  held  constant  at  =  2.4  E/min. 

Instrumentation  signals  generated  on  the  table  were  acquired  through  a  custom- 
built  forty-channel  instrumentation  slip  ring  using  a  data  acquisition  system. 
Temperatures,  mass  flow  rates,  accelerations,  and  voltages  were  all  measured  using  a 
data  acquisition  mainframe  (Agilent  VXI  E8408A)  with  a  command  module  (Agilent 
E1406A),  514  digit  multimeter  module  (Agilent  E1411B),  and  a  64-channel  3-wire 
multiplexer  module  (Agilent  E1476A).  The  rotational  speed  of  the  centrifuge  table, 
heater  power,  and  other  low  voltage  control  devices  on  the  table  were  controlled  using  an 
8/1 6-channel  D/A  converter  module  (Agilent  E1418A).  Communication  between  the  data 
acquisition  unit  and  the  computer  was  established  using  a  general  purpose  interface  bus 
(GPIB)  coupled  with  a  custom-designed  EabVIEW  virtual  instrument. 
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Gathering  temperature  data  from  rotating  machinery  using  slip  rings  presents 
unique  problems.  First,  when  the  thermocouple  wires  are  connected  to  the  wires  leading 
to  a  slip  ring,  at  least  one  extra  junction  is  created,  depending  on  the  materials  of  the 
thermocouple  wires.  To  avoid  this  problem,  a  Type  E  thermocouple  amplifier  was 
installed  on  the  centrifuge  table  (Omega  OM7-47-E-07-2-C)  with  internal  cold  junction 
compensation.  This  converted  the  millivoltage  signals  from  the  thermocouples  to  0  to  10 
V  signals  without  the  creation  of  extra  junctions.  Another  problem  that  is  present  when 
slip  rings  are  used  is  electrical  noise.  This  problem  was  reduced  (not  eliminated)  by  the 
use  of  a  low-pass  filter  for  each  of  the  thermocouple  signals  coming  from  the  table  before 
the  data  acquisition  system. 

The  test  article,  a  titanium-water  loop  heat  pipe,  was  developed  for  AERL/RZPS 
by  Advanced  Cooling  Technologies  (ACT),  Inc.,  in  Eancaster,  PA,  under  contract 
PA8601-06-P-0076.  Initial  design  parameters  set  by  AERL/RZPS  were  to  develop  a  loop 
heat  pipe  capable  of  a  minimum  heat  load  of  500  W  and  minimum  heat  flux  of  3  W/cm  . 
The  minimum  transport  line  length  was  2  m  to  simulate  relevant  aircraft  geometries.  An 
evaporator  operating  temperature  of  200°C  and  condenser  operating  temperature  between 
5  and  140°C  were  selected  to  match  relevant  acquisition  and  rejection  temperatures 
aboard  aircraft.  The  evaporator  and  condenser  dimensions  were  selected  to  be  20.32  x 
10.16  cm  and  30.48  x  28.56  cm,  respectively,  to  match  commercial  off-the-shelf  heaters 
and  cold  plates.  A  summary  of  the  requested  design  parameters  can  be  seen  in  Table  2.1. 
After  several  design  iterations,  ACT  delivered  the  loop  heat  pipe  shown  in  Eigure  2.4.  A 
summary  of  the  loop  heat  pipe  specifications  can  be  seen  in  Table  2.2.  The  LHP  was 
instrumented  with  twelve  type  E  exposed  tip  thermocouples  as  seen  in  Eigure  2.5.  A 
summary  of  their  locations  can  be  seen  in  Table  2.3. 

The  loop  heat  pipe  was  mounted  onto  the  centrifuge  table  such  that  the  centerline 
of  the  tubing  coincided  with  the  outer  table  radius  as  much  as  possible.  Small  deviations 
existed  since  the  condenser  section  and  the  evaporator/compensation  chamber  were  both 
straight.  This  induced  a  non-uniform  radial  acceleration  field  over  the  lengths  of  these 
sections  that  needed  to  be  quantified.  Stands  were  designed  using  G-7  phenolic  to  mount 
the  loop  heat  pipe  with  support  at  the  compensation  chamber,  evaporator,  condenser,  and 
transport  lines  (Eigure  2.6).  The  tops  of  these  stands  were  anchored  to  the  table  to  reduce 
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deflection  when  the  table  was  rotating.  A  survey  was  taken  at  22  locations  on  the  loop 
heat  pipe  to  determine  how  far  various  portions  of  the  loop  heat  pipe  were  from  the 
centerline  radius.  The  loop  heat  pipe  had  a  minimum  radius  to  centerline  of  119.2  cm 
and  a  maximum  radius  to  centerline  of  123.3  cm.  The  entire  loop  heat  pipe  fitted  within 
4.6  cm  for  a  percent  acceleration  difference  of  3.8%.  Complete  survey  data  can  be  seen 
in  Appendix  D.  To  minimize  heat  loss  to  the  environment,  the  entire  assembly  was 
thoroughly  insulated  using  Kaowool  blankets  and  aluminum  foil.  The  assembly  was 
placed  inside  an  aluminum  frame  (80/20,  Inc.)  for  structural  support  and  enclosed  with 
sheet  metal  sides  to  minimize  convective  heat  losses. 

During  operation,  heat  was  applied  to  the  LHP  at  the  evaporator  while  the  heat 
transfer  to  the  compensation  chamber  was  independently  controlled.  A  mica  heater 
(Minco)  was  located  between  the  evaporator  body  and  a  ceramic  fiber  insulative  layer, 
followed  by  the  evaporator  stand.  A  flexible  electric  heat  tape  (Thermolyne)  was  wound 
around  the  compensation  chamber  and  surrounded  by  Kaowool  insulation  and  aluminum 
foil  to  minimize  heat  losses.  In  normal  operation,  the  compensation  chamber  is  not 
insulated  and  the  temperature  is  closely  controlled  during  operation.  For  these 
experiments,  insulating  the  LHP,  including  the  compensation  chamber,  was  selected  to 
mimic  a  typical  configuration  of  a  LHP  in  an  aircraft  environment  where  bay 
temperatures  could  be  higher  than  the  LHP  temperatures.  This  would  minimize  parasitic 
heat  gain,  and  reduce  the  use  of  external  heaters  or  coolers  on  the  compensation  chamber. 
As  a  result,  the  LHP  compensation  chamber  was  allowed  to  “float”  into  equilibrium  with 
the  evaporator  and  condenser,  rather  than  controlling  the  temperature  of  the  evaporator 
by  controlling  the  compensation  chamber  temperature. 

As  previously  mentioned,  the  centrifuge  table  was  equipped  with  an  on-board 
fluid  loop  for  dissipating  heat  from  sources  on  the  table,  which  used  ethylene  glycol  as  its 
working  fluid.  In  the  present  experiment,  it  was  desired  to  have  the  option  of  operating 
the  LHP  condenser  section  at  elevated  temperatures,  so  a  high-temperature  fluid  loop  was 
constructed  and  mounted  to  the  centrifuge  table  to  act  as  an  interface  between  the  LHP 
and  the  low-temperature  fluid  loop,  as  shown  in  Figure  2.7.  The  high-temperature 
working  fluid  (Brayco  Micronic  889  polyalphaolefin  or  PAO  oil)  flowed  from  the 
custom-made  copper  reservoir  into  a  positive  displacement  gear  pump  (Tuthill).  After 
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passing  through  a  filter  and  a  flow-straightening  section,  the  PAO  was  directed  through 
the  turbine  flow  meter  (Omega  FTB-9506).  An  electrical  tape  heater  was  mounted  to  the 
copper  tubing  after  the  flow  meter  to  allow  for  preheating  the  PAO  prior  to  reaching  the 
calorimeter  on  the  condenser  section,  which  consisted  of  three  heat  exchangers  plumbed 
in  series  and  mounted  to  the  condenser  section.  Type  E  thermocouple  probes  were 
installed  at  the  inlet  and  outlet  of  the  three  heat  exchangers  for  calorimetry  (TCOO  and 
TCOl),  and  another  was  placed  prior  to  the  flow  meter  (TC03).  This  was  needed  due  to 
the  dependence  of  the  viscosity  of  PAO  on  temperature.  After  the  PAO  exited  the  three 
heat  exchangers  on  the  condenser,  it  flowed  to  a  liquid/liquid  heat  exchanger  that 
transferred  heat  from  the  high-temperature  coolant  loop  to  the  low-temperature  ethylene 
glycol  loop.  The  PAO  then  returned  to  the  reservoir. 

Four  grounded  probe  thermocouples  for  the  high  temperature  loop  and  twelve 
exposed  tip  type  E  thermocouples  mounted  on  the  EHP  were  used  in  the  experiment. 
Thermocouple  calibrations  were  conducted  over  two  temperature  ranges  depending  on 
the  anticipated  operating  temperatures.  The  grounded  probe  thermocouples  were  used  for 
calorimetry,  coolant  flow  meter  calibration  and  the  measurement  of  the  ambient 
temperature,  where  the  error  needed  to  be  minimized.  These  four  thermocouples  were 
calibrated  over  the  anticipated  range  of  20  to  145°C  in  5°C  intervals.  The  twelve 
exposed  tip  thermocouples  were  mounted  on  the  EHP  in  various  locations  and  needed  to 
be  calibrated  over  the  full  range  of  20  to  230°C  in  5°C  increments.  The  calibration 
procedure  consisted  of  using  two  separate  recirculating  chiller  baths  (Brinkmann  Lauda 
RCS  20-D,  r  =  20  to  140°C;  Hart  Scientific  6330,  T  =  40  to  230°C)  with  PAO  as  the 
working  fluid  to  achieve  the  required  temperature  range.  The  temperature  readings  from 
the  sixteen  thermocouples  were  compared  to  a  NIST-traceable  platinum  resistance 
temperature  detector  (Hart  Scientific  RTD  1502A)  with  a  resolution  of  ±  0.009°C.  To 
ensure  that  the  bath  had  reached  steady  state  at  a  given  temperature,  the  RTD  temperature 
was  continuously  monitored.  When  the  standard  deviation  of  100  readings  dropped 
below  the  specified  threshold  of  0.005°C,  100  readings  from  the  thermocouples  were 
sampled,  stored  in  an  array,  and  the  bath  temperature  was  changed.  For  repeatability,  the 
bath  temperature  was  first  incremented  from  the  lowest  temperature  to  the  highest 
temperature,  and  then  decremented  from  highest  to  lowest,  and  the  two  sets  of  100  data 
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points  collected  for  each  thermocouple  at  a  given  temperature  were  used  to  determine 
two  average  readings.  Plots  of  the  RTD  temperature  versus  each  thermocouple 
temperature  were  generated,  and  polynomial  trend  lines  were  fitted  for  each 
thermocouple  as  can  be  seen  in  Appendix  C.  A  fifth-order  polynomial  was  selected  since 
it  reduced  the  maximum  deviation  from  the  data  by  approximately  a  factor  of  four  over  a 
first-order  trend  line.  The  uncertainty  associated  with  each  thermocouple  was  determined 
by  accounting  for  four  sources  of  error:  the  stated  uncertainty  of  the  RTD,  the  confidence 
interval  of  the  RTD  average  reading  at  a  confidence  level  of  0.95,  the  confidence  interval 
of  the  thermocouple  average  reading  at  a  confidence  level  of  0.95,  and  the  maximum 
deviation  of  the  temperature  calculated  using  the  polynomial  curve  fit  from  the  actual 
measured  temperature. 

The  turbine  flow  meter  used  in  the  high-temperature  fluid  loop  was  calibrated  to 
achieve  accurate  results  for  the  amount  of  heat  extracted  from  the  LHP.  This  was  critical 
for  the  calculation  of  the  evaporative  heat  transfer  coefficient  and  the  thermal  resistance 
of  the  LHP.  Since  the  viscosity  of  the  PAO,  used  in  the  high-temperature  fluid  loop, 
changes  significantly  with  temperature,  a  “calibration  surface”  was  generated  that  related 
the  output  voltage  of  the  flow  meter  and  the  temperature  of  the  PAO  at  the  entrance  of 
the  flow  meter  to  the  mass  flow  rate.  The  calibration  setup  consisted  of  a  recirculating 
chiller  bath  (Brinkmann  Lauda  RCS  20-D)  filled  with  PAO  from  the  same  source  as  used 
in  the  high-temperature  fluid  loop.  The  gear  pump,  inline  filter  (Whitey  SS56S6  140 
micron)  and  a  calibrated  grounded  thermocouple  probe,  from  the  high-temperature  fluid 
loop,  were  installed  in  a  line  from  the  bath  to  the  turbine  flow  meter  (Omega  FTB-9506) 
and  signal  conditioner  (Omega  FLSC-61).  Flow  straightening  sections  upstream  and 
downstream  were  placed  according  to  the  manufacturer's  instructions.  A  three-way  valve 
was  installed  after  the  flow  meter,  which  allowed  the  entire  flow  system  to  reach  a  steady 
temperature.  Once  the  temperature  was  steady,  the  flow  was  diverted  to  a  catch  basin  for 
a  specified  amount  of  time.  The  voltage  from  the  flow  meter  and  the  temperature  from 
the  thermocouple  were  recorded  during  this  time,  and  when  the  basin  was  full,  the  flow 
was  again  diverted  to  recirculating  the  PAO  back  to  the  chiller  bath.  All  of  the  data  was 
collected  through  the  instrumentation  slip  rings  on  the  centrifuge  table  to  the  data 
acquisition  system  to  capture  all  errors  inherent  to  the  centrifuge  table  test  bed.  A  lab 
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scale  (Mettler  PC4400)  was  used  to  determine  the  mass  eolleeted  during  a  given  test  run 
to  within  ±0.3  gm.  During  eaeh  measurement,  as  many  data  points  as  possible  were 
eolleeted  aeross  the  time  span  with  the  limiting  faetor  being  the  iteration  time  on  the 
Lab  VIEW  software.  The  minimum  number  of  data  points  eolleeted  for  any  given  run 
was  437.  The  voltages  and  temperatures  were  averaged  and  a  eonfidenee  interval  was 
ealeulated  based  on  a  eonfidenee  level  of  0.95  for  eaeh  test  run.  The  test  was  repeated  for 
a  total  of  five  averaged  data  points  for  eaeh  nominal  temperature  and  flow  rate.  These 
tests  were  eompleted  over  the  range  of  T  =  20  to  120°C  in  intervals  of  25°C  and  flow 
rates  ranging  from  m  =  0.0064  to  0.025  kg/s  in  intervals  of  approximately  0.002  kg/s.  A 
3-D  paraboloid  regression  equation  was  generated  using  SigmaPlot  to  relate  temperature, 
flow  meter  voltage,  and  mass  flow  rate,  and  was  given  by 

rhcp  =  yo  +  aT  +  bV  +  cT^  +  dV^  (2.1) 

where  yo,  a,  b,  c,  and  d  are  ealibration  eonstants  (Appendix  C).  The  general  root-sum- 
square  uneertainty  equation  used  for  all  uneertainties  was  given  by 

where  y  =J{x\,  xi,  ...).  The  uneertainty  of  the  mass  flow  rate  measurement  was  affeeted 
by  the  maximum  deviation  of  the  regression  equation  from  the  aetual  data,  the  eonfidenee 
interval  for  the  temperature  and  flow  meter  voltage  measurements,  the  root-sum-square 
total  error  assoeiated  with  the  seale  and  stopwateh  given  by 


and  the  root-sum-square  error  assoeiated  with  the  temperature  and  voltage  measurements 
given  by 

r  2  21^/^ 

Arhv/T  =  [((a  +  2cT)AT)  +  ((h  +  2dV)AV)  J  (2.4) 
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The  percent  error  on  the  mass  flow  rate  decreased  with  increasing  flow  rate.  Since  the 
mass  flow  rate  was  kept  constant  at  =  0.0077  kg/s,  the  uncertainty  associated  with 

that  setting  was  4.0%. 

The  heat  transferred  from  the  LHP  condenser  to  the  cold  plate,  gout,  was  defined 
as 


l3out  ~  ^cp^p,PAo(^out  ^in)  (2.5) 

A  linear  fit  equation  for  Cpp^p  as  a  function  of  temperature  was  developed  by  Ghajar  et 

al.  (1994)  and  used  in  equation  (2.5)  (Appendix  E).  The  average  evaporative  heat 
transfer  coefficient  was  defined  as 


h  = 


Qout 

nDLif,  -  T^) 


(2.6) 


where  D  is  the  inside  diameter  of  the  evaporator  shell,  L  is  the  length  of  the  evaporator, 
Eg  is  the  average  evaporator  temperature  measured  by  the  four  thermocouples  embedded 
in  the  wall  between  the  heater  and  the  wick  (Figure  2.5(b)),  and  TV  is  the  external 
temperature  of  the  vapor  line  at  the  outlet  of  the  evaporator.  The  heat  rejected  to  the  cold 
plate,  gout,  was  selected  as  it  was  the  best  estimate  of  heat  actually  transported  by  the 
LHP.  The  thermal  resistance  of  the  loop  heat  pipe,  R,  was  determined  using  the  average 
evaporator  temperature  and  the  average  temperature  of  the  cold  plate,  and  was  defined  as 


R  = 


T  —T 

^  e  ^  cp 
l2out 


(2.7) 


where  is  the  average  cold  plate  temperature.  The  root-sum-square  uncertainty  of  gout, 
h  ,  and  R  are  given  by 


^Qout  —  [(6'p,pAo(2’out  2’jn)Arh(.p)  +  (fiicpC^’out  2’j„)ACppAo) 

2  2  ^ 
(^cp6'p,PA0^2)m()  +  (~jficp6'p,PA0^2jf,)  ] 
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The  uncertainty  of  Cpp^o  was  estimated  by  Ghajar  et  al.  to  be  0.5%  of  the  value.  For 

each  steady  state  condition,  151  data  points  were  collected  from  each  sensing  device 
representing  five  minutes  of  data.  Measured  values  were  averaged  and  uncertainties  were 
calculated  based  on  the  fixed  error  of  each  instrument  and  the  confidence  interval  for  the 
average  at  a  confidence  level  of  0.95.  A  summary  of  the  uncertainties  for  this  experiment 
can  be  found  in  Table  2.4.  Details  of  the  uncertainty  analysis  can  be  found  in  Appendix 
B. 

2.4.  Results  and  Discussion 

The  purpose  of  this  series  of  experiments  was  to  determine  the  operating 
characteristics  of  a  titanium-water  loop  heat  pipe  subjected  to  changes  in  evaporator  heat 
input,  compensation  chamber  heat  input,  and  radial  acceleration.  Steady  state  and 
transient  temperature  data  were  collected  which  provided  insight  into  the  fluid-thermal 
behavior  of  the  LHP.  The  raw  data  was  reduced  to  obtain  the  evaporative  heat  transfer 
coefficient,  thermal  resistance,  and  evaporator  wall  superheat  in  terms  of  the  heat 
transported  and  radial  acceleration  level.  Quasi-steady  phenomena  and  dry-out  of  the 
LHP  were  observed  and  quantified  in  a  performance  map. 

Figure  2.8  presents  a  typical  stationary  (uz  =  1.0  g,  Or  =  0.0  g)  cold-start  test  of  the 
LHP,  which  consisted  of  the  following:  With  the  LHP  at  ambient  conditions,  the 
recirculating  chiller  in  the  low-temperature  fluid  loop  was  set  to  Teg  =  35°C.  Heat  was 
applied  as  a  step  function  to  the  evaporator  section  (in  this  case,  Qin  =  600  W)  while  the 
pump  for  the  high-temperature  fluid  loop  was  simultaneously  turned  on  (m^^  =  0.0077 
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kg/s).  The  mass  flow  rate  of  the  high-temperature  fluid  loop  was  maintained  constant  at 
this  value  throughout  this  series  of  experiments  to  minimize  the  uncertainty  associated 
with  the  calorimetry  of  the  cold  plate.  Figure  2.8(a)  shows  the  transient  temperature 
response  of  the  evaporator,  vapor  line,  and  calorimeter  inlet  and  outlet.  The  temperatures 
appear  to  become  steady  after  approximately  6000  s.  However,  in  order  to  determine 
when  steady  state  occurred  the  time  rate  of  change  of  the  temperatures  was  averaged  over 
15  min,  intervals  and  plotted  with  respect  to  time  as  shown  in  Figure  2.8(b).  It  was 
observed  that  dT/dt  approached  zero  shortly  after  6000  s,  but  for  times  greater  than  6000 
s,  significant  oscillations  occurred.  The  oscillations  in  ATI  At  were  not  apparent  in  the  raw 
temperature  traces,  but  steady  state  was  found  to  occur  at  approximately  18,000  s.  This 
was  further  demonstrated  by  calculating  the  thermal  resistance  and  heat  transfer 
coefficient  for  this  test  at  different  times,  as  shown  in  Figure  2.8(c).  This  methodology 
was  used  throughout  testing  to  ensure  that  a  repeatable  steady  state  was  reached. 

Figure  2.9  also  shows  transient  temperature  traces  during  the  Qin  =  600  W  test 
described  in  the  previous  paragraph.  In  Figure  2.9(a),  the  evaporator  temperature 
(TC04)  increased  very  quickly  while  the  rest  of  the  LHP  did  not  react.  After 
approximately  60  s,  the  thermocouple  located  on  the  vapor  line  nearest  to  the  exit  of  the 
evaporator  (TC08)  suddenly  increased.  This  was  followed  in  turn  by  increases  in 
temperature  reflected  by  the  thermocouples  located  throughout  the  condenser  section. 
This  shows  the  progression  of  the  saturated  vapor  clearing  the  condenser  section  of 
liquid,  which  was  subsequently  displaced  into  the  evaporator  section  and  the 
compensation  chamber  via  the  bayonet  tube.  Figure  2.9(b)  shows  that  the  evaporator 
temperature  was  significantly  higher  than  the  condenser  temperatures,  which  led  to  a 
relatively  high  value  of  thermal  resistance,  which  will  be  discussed  in  detail  below. 

Figure  2.10  shows  temperature  traces  in  the  condenser  (TC09  through  TCI 3)  and 
at  the  bayonet  inlet  (TC14).  Each  figure  shows  the  transient  temperature  after  the 
stationary  LHP  reached  steady  state  conditions  at  heat  inputs  ranging  from  100  ^  Qm  ^ 
600  W.  In  Figure  2.10(a),  with  =  100  W,  the  liquid  entering  the  bayonet  tube  was 
highly  subcooled  at  approximately  40°C.  At  this  heat  input  level,  the  majority  of  the 
condenser  was  flooded  with  subcooled  liquid.  In  fact,  only  TC09  (condenser  inlet) 
indicated  two-phase  flow.  Figure  2.10(b),  with  gin  =  200  W,  was  a  unique  case  that  is 
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described  further  in  the  following  paragraph  (Figure  2.11).  Figure  2.10(c)  to  Figure 
2.10(f)  shows  that  the  2(p-l(p  point  progressed  through  the  liquid  line  as  heat  input 
increased  until  it  reached  the  bayonet  inlet.  If  the  heat  input  at  the  evaporator  is  high 
enough,  saturated  vapor  will  pass  through  the  bayonet  tube  and  reach  the  evaporator 
section.  This  point  represents  a  performance  limit  to  the  LHP  operation  because  if  vapor 
enters  the  evaporator,  the  wick  will  dry  out  and  the  LHP  will  overheat. 

Figure  2.1 1  shows  the  oscillatory  behavior  of  the  LHP  for  the  heat  input  of  Q\a  = 
200  W.  Initially,  at  t  =  0,  the  evaporator  temperatures  (TC04,  TC05,  TC06,  and  TC07) 
ranged  from  66  to  68°C.  The  evaporator  temperature  nearest  to  the  bayonet  tube  outlet 
(TC07)  was  the  lowest,  which  indicated  that  the  subcooled  liquid  that  entered  the 
evaporator  tended  to  reduce  the  evaporator  temperature  at  this  point.  The  vapor  line  and 
condenser  temperatures  (TC08  through  TC13)  ranged  from  46  to  58°C.  The  vapor  line 
(TC08)  was  the  highest,  with  the  first  three  thermocouples  in  the  condenser  (TC09, 
TCIO,  TCI  1)  decreasing  slightly.  The  vapor  became  saturated  within  the  condenser,  and 
condensation  formed  on  the  interior  walls  of  the  tubing.  From  the  point  at  which  the 
quality  of  the  working  fluid  was  x  =  1  (saturated  vapor)  to  where  it  reached  x  =  0 
(saturated  liquid),  the  temperature  should  have  been  constant,  except  for  the  fact  that  the 
pressure  dropped  slightly  due  to  viscous  losses.  This  drop  in  the  saturation  pressure  in 
turn  decreased  the  saturation  temperature.  Past  TCll,  the  other  condenser  temperatures 
(TC12,  TC13)  dropped  significantly.  This  showed  that  the  2(p-lq)  point,  where  x  =  0, 
occurred  between  TCll  and  TC12.  The  working  fluid  after  this  point  became  a 
subcooled  liquid,  where  the  temperature  drop  was  due  to  sensible  heat  extraction  by  the 
cold  plates.  Interestingly,  at  t  =  0,  the  temperature  at  the  bayonet  inlet  (TC14)  was  higher 
than  the  outlet  of  the  condenser.  Under  typical  operation,  this  was  not  the  case  due  to 
convective  losses  from  the  liquid  lines. 

As  time  progressed  from  t  =  0  (Figure  2.11(a)),  several  things  occurred  nearly 
simultaneously.  The  evaporator  thermocouple  nearest  to  the  vapor  manifold  (TC07) 
suddenly  decreased,  which  indicated  movement  of  subcooled  liquid  from  the  exit  of  the 
bayonet  tube  into  the  evaporator.  The  junction  between  the  evaporator  and  the 
compensation  chamber  (TCI 5)  increased  and  then  decreased  in  temperature  over  a 
relatively  short  period.  This  was  due  to  warm  liquid  in  the  evaporator  section  being 
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pushed  through  the  grooves  into  the  compensation  chamber,  followed  by  cooler  liquid 
from  the  bayonet  tube  exit.  The  inlet  of  the  bayonet  tube  (TCI 4)  decreased,  and  the  two 
thermocouples  measuring  the  subcooled  liquid  in  the  condenser  increased  (TC12  and 
TCI  3).  Again,  this  was  indicative  of  movement  of  the  slug  of  liquid  that  existed  from  the 
2q)-l(p  point  in  the  condenser  to  the  meniscus  within  the  grooves  of  the  secondary  wick 
inside  the  evaporator  section,  as  shown  in  Figure  2.2(a).  The  dramatic  increase  in  the 
condenser  section  (TC12)  shows  that  the  2(p-l(p  point  moved  from  between  TCll  and 
TC12,  across  the  TC12  location,  and  then  between  TC12  and  TC13  as  shown 
schematically  in  Figure  2.11(b).  In  fact,  TCI 2  increased  to  the  saturated  vapor 
temperature  existing  within  the  first  half  of  the  condenser. 

At  approximately  t  =  80  s,  the  temperatures  in  the  evaporator  and  the  bayonet  tube 
inlet  (TC14)  started  to  increase,  while  the  condenser  temperatures  TC12  and  TC13 
decreased.  This  behavior  indicated  that  the  liquid  slug  had  reversed  direction;  i.e.  the  2(p- 
l(p  point  re-crossed  thermocouple  location  TC12  in  the  condenser.  The  significant  rise  in 
the  bayonet  inlet  temperature  TCI 4  shows  that  warm  liquid  originally  in  the  evaporator 
was  now  flooding  back  through  the  bayonet  tube  into  the  liquid  line.  This  movement  of 
liquid  out  of  the  evaporator  may  be  due  to  the  sudden  appearance  of  a  vapor  bubble 
within  the  wick  structure  of  the  evaporator  section  which  would  tend  to  drive  the  heated 
liquid  in  the  evaporator  in  the  opposite  direction.  As  can  be  seen  in  Figure  2.11(a),  the 
period  of  the  oscillation  was  approximately  150  s.  This  type  of  percolation  is  not  typical 
of  a  fully  operational  LHP,  but  is  actually  closer  to  the  behavior  of  a  pulsating  heat  pipe. 
Discussion  of  flow  reversal  within  LHPs  in  the  literature  was  limited  to  startup  and 
shutdown  operation.  Douglas  et  al.  (1999)  discussed  flow  reversal  in  LHPs  as  a 
phenomenon  that  occurred  during  startup  and  continued  until  the  capillary  pressure  in  the 
secondary  wick  could  no  longer  maintain  the  system  pressure  drop.  Cimbala  et  al.  (2004) 
used  neutron  radiography  to  visualize  LHP  operation  and  observed  flow  reversal  only 
occurred  when  the  heat  input  was  reduced  to  =  0  W.  It  was  concluded  that  with  no 
heat  input,  convective  and  radiative  heat  transfer  from  the  LHP  to  the  ambient  caused  the 
flow  reversal.  In  general,  flow  reversal  was  not  discussed  as  part  of  normal  operation. 
However,  in  the  present  experiment,  flow  reversal  was  found  at  some  operating  points. 
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Figure  2.12  shows  the  various  steady  state  LHP  temperatures  versus  transported 
heat  for  the  stationary  ease.  The  four  evaporator  temperatures  in  Figure  2.12(a)  inereased 
monotonieally  with  heat  transported,  but  diverged  from  the  vapor  outlet  temperature. 
The  behavior  of  the  condenser  temperatures  with  heat  transported  was  slightly  different, 
as  shown  in  Figure  2.12(b).  At  the  lowest  heat  input  value  (gm  =  100  W),  a  significant 
temperature  drop  was  present  between  the  inlet  of  the  condenser  (TC09)  and  the 
thermocouples  within  the  condenser.  This  shows  that  the  2q)-l(p  point  resided  between 
TC09  and  TCIO,  which  means  that  very  little  of  the  available  condenser  was  being  used 
for  two-phase  condensation.  This  condition  also  shows  that  the  liquid  returning  to  the 
evaporator  section  (TC14)  was  highly  subcooled.  As  the  heat  input  increased  to  Qin  = 
200  W,  the  temperatures  measured  at  TCIO  and  TCI  1  rose  to  match  that  at  TC09,  which 
means  that  the  time  averaged  location  of  the  2(p-l(p  point  moved  farther  into  the 
condenser  (between  TCll  and  TC12).  At  a  heat  input  of  gm  =  300  W,  the  2(p-l(p  point 
traveled  past  the  end  of  the  condenser  into  the  liquid  lines  such  that  all  of  the  condenser 
temperatures  matched  the  evaporator  outlet  temperature  (TC08).  As  the  heat  input 
increased,  the  condenser  temperatures  continued  to  rise.  However,  the  evaporator  outlet 
temperature  increased  at  a  faster  rate,  which  is  indicative  of  an  increased  superheat 
penalty. 

Figure  2.13  shows  the  thermal  performance  of  the  stationary  LHP  for  heat  inputs 
ranging  from  Qi^  =  100  to  600  W.  The  evaporative  heat  transfer  coefficient.  Figure 
2.13(a),  decreased  monotonieally  with  transported  heat.  This  behavior  was  controlled  by 
the  slope  of  the  average  evaporator  temperature  versus  that  of  the  evaporator  outlet,  as 
shown  in  Figure  2.12(a).  The  temperature  difference  -T^')  defined  in  equation  (2.6) 

increased  more  rapidly  than  Qout,  which  resulted  in  an  overall  decrease  in  h  .  As  dry-out 
was  approached,  more  of  the  wick  in  the  evaporator  section  was  depleted  of  liquid,  which 
tended  to  increase  the  evaporator  temperature.  The  thermal  resistance  of  the  stationary 
LHP  versus  heat  transported  is  presented  in  Figure  2.13(b),  where  it  is  seen  to  decrease, 
reach  a  minimum,  and  then  increase.  At  low  power  inputs,  the  relatively  large 
temperature  drop  defined  by  equation  (2.7),  ~Tcp),  drives  the  thermal  resistance  to  a 
high  value.  This  temperature  drop  was  a  result  of  the  fact  that  most  of  the  condenser 
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section  was  flooded  by  subcooled  liquid  which  was  close  to  the  cold  plate  temperature. 
As  the  2(p-lq)  point  moved  through  and  then  exited  the  condenser,  the  temperature  drop 
decreased  with  transported  heat,  which  decreased  the  thermal  resistance.  The  minimum 
R  corresponds  to  the  point  in  Figure  2.12(b)  where  the  2(p-l(p  point  just  exited  the 
condenser.  Past  this  point,  the  evaporator  section  increased  in  temperature  more  rapidly 
than  the  condenser  section,  which  resulted  in  the  thermal  resistance  increasing  with 
transported  heat.  The  wall  superheat,  defined  as  the  difference  between  the  average 
evaporator  temperature  and  the  temperature  of  the  evaporator/compensation  chamber 
junction,  was  found  to  monotonically  increase  with  an  increasing  amount  of  transported 
heat.  With  respect  to  the  evaporative  heat  transfer  coefficient,  thermal  resistance,  and 
wall  superheat,  no  notable  difference  was  observed  between  starting  the  LHP  while  the 
unit  was  at  ambient  temperature  versus  a  step  change  in  the  evaporator  heat  input  from  a 
lower  to  higher  value  or  a  higher  to  lower  value.  A  summary  of  the  stationary  steady 
state  data  points  and  the  path  to  reach  steady  state  can  be  seen  in  Table  2.5. 

Figure  2.14  shows  the  operating  characteristics  and  performance  of  the  stationary 
LHP  for  an  evaporator  heat  input  of  Qi^  =  500  W  while  varying  the  compensation 
chamber  heat  input  from  =  0  to  50  W.  For  this  particular  test,  the  LHP  was  allowed 
to  achieve  steady  state  conditions  for  the  given  evaporator  heat  input,  after  which  the 
compensation  chamber  heat  input  was  incremented  in  steps  of  5  W.  In  Figure  2.14(a) 
and  Figure  2.14(b),  for  gcc  =  0  W,  the  evaporator  temperatures  were  relatively  uniform, 
where  the  vapor  exiting  the  evaporator  was  slightly  superheated  and  the  2(p-l(p  point  was 
out  of  the  condenser.  When  a  small  amount  of  heat  was  input  to  the  system  through  the 
compensation  chamber  (gcc  =  5  W),  the  evaporator  temperatures  and  the  evaporator  exit 
temperature  both  decreased  while  the  condenser  temperatures  remained  constant.  This 
trend  continued  until  approximately  gcc  =  15  W,  at  which  point  the  evaporator 
temperature  leveled  off,  the  evaporator  exit  temperature  decreased  to  the  saturation 
temperature  within  the  condenser,  and  the  condenser  outlet  temperature  dropped  below 
the  saturation  temperature.  The  decrease  in  the  average  evaporator  temperature 
significantly  affected  the  evaporative  heat  transfer  coefficient  and  the  thermal  resistance, 
as  shown  in  Figure  2.14(c).  In  fact,  h  increased  by  68%  with  an  increase  in  the  overall 
heat  input  of  only  3%.  The  drop  in  the  condenser  outlet  temperature  indicated  that  the 
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2(p-l(p  point  moved  from  the  liquid  line  into  the  eondenser  seetion.  For  gcc  >  20  W,  the 
evaporator  temperatures  inereased  and  the  2(p-l(p  point  eontinued  to  move  toward  the 
evaporator  whieh  resulted  in  an  inerease  in  the  thermal  resistanee.  Ku  (1999)  indieated 
that  operating  the  eompensation  ehamber  at  a  higher  temperature  by  using  an  external 
heater  in  effeet  inereases  the  amount  of  subeooling  in  the  eondenser  and  liquid  return 
line.  Aeeording  to  Ku  (1999),  this  subeooling  is  neeessary  to  balanee  the  additional  heat 
input  and  results  in  underutilizing  the  eondenser  and  a  degradation  of  the  thermal 
eonduetanee.  In  the  present  experiment,  this  eonelusion  held  true  for  Q^c  >  20  W,  as  the 
amount  of  subeooling  inereased  the  thermal  resistanee  and  deereased  the  evaporative  heat 
transfer  eoeffieient  by  way  of  an  inereased  superheat  penalty.  This  did  not  hold  true  for 
Qcc  <  20  W.  When  the  LHP  operated  at  Q\n  =  500  W,  the  2(p-l(p  interfaee  was  loeated  in 
the  liquid  return  line.  Inereasing  the  heat  input  to  the  eompensation  ehamber  moved  the 
2(p-l(p  interfaee  to  the  eondenser  outlet  at  Qcc  =  15  W.  Operation  at  this  point  maximized 
the  amount  of  heat  transfer  due  to  eondensation  with  the  added  benefit  of  eooler  liquid  in 
the  eompensation  ehamber  and  evaporator  whieh  deereased  the  thermal  resistanee  and 
inereased  the  evaporative  heat  transfer  eoeffieient. 

Also  of  interest  is  the  temperature  inerease  at  the  bayonet  inlet  (TCI 4)  starting 
when  Qcc  =  35  W  seen  in  Figure  2.14(b).  Figure  2.15  shows  the  transient  temperature 
traees  of  the  eondenser,  bayonet  tube,  and  evaporator/eompensation  ehamber  junetion  for 
Qcc  =  25  to  50  W.  In  Figure  2.15(a)  and  Figure  2.15(b),  with  Qcc  =  25  and  30  W, 
subeooled  liquid  moved  through  the  bayonet  inlet  as  seen  in  typieal  operation.  In  Figure 
2.15(e)  through  Figure  2.15(f),  with  Qcc  =  35  through  50  W,  a  sudden  inerease  in 
temperature  at  the  bayonet  inlet  (TCI 4)  showed  that  flow  reversal  oeeurred  in  the 
evaporator  seetion.  This  was  similar  to  the  oseillating  phenomena  deseribed  for  Qi^  = 
200  W  exeept  that  the  liquid-vapor  meniseus  in  the  seeondary  wiek  was  driven  baekward 
by  the  elevated  vapor  pressure  within  the  eompensation  ehamber,  whieh  was  due  to  the 
heat  input  at  the  shell  of  the  eompensation  ehamber.  In  addition,  the  temperature  of  the 
evaporator/eompensation  ehamber  interfaee  (TCI 5)  did  not  vary  appreeiably,  whieh  was 
different  than  that  seen  at  gin  =  200  W.  The  liquid-vapor  meniseus  moved  baekward  due 
to  the  inereased  pressure  within  the  eompensation  ehamber  until  a  point  at  whieh  the 
pressure  was  balaneed.  Forward  flow  then  resumed  and  heat  was  lost  through  the  liquid 
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line,  shown  by  the  slow  deerease  in  temperature  at  the  bayonet  inlet  (TC14)  to  the 
temperature  of  the  subeooled  liquid  at  the  condenser  outlet  (TCI 3).  A  summary  of  the 
steady  state  data  points  for  gcc  =  0  to  50  W  can  be  seen  in  Table  2.6. 

To  further  explore  compensation  chamber  heat  input  and  heat  loss  to  the  ambient, 
Table  2.7  shows  the  effect  of  operating  the  stationary  LHP  for  Qi^  =  500  W  with  the 
compensation  chamber  uninsulated,  insulated,  temperature  controlled  to  Tec  =  72.8°C  via 
simultaneous  heat  input  to  the  compensation  chamber  (Qcc  =  20  W)  and  evaporator,  and 
preconditioning  the  temperature  to  Tec  =  72.8°C  via  heat  input  (Qcc  =  100  W  reduced  to 
Qcc  =  20  W)  prior  to  heat  input  to  the  evaporator.  For  this  series  of  tests,  thermocouple 
TCI 5  was  relocated  to  the  top  side  of  the  compensation  chamber  to  directly  monitor  its 
operating  temperature.  It  was  observed  that  the  average  evaporator  temperature 
increased,  the  evaporative  heat  transfer  coefficient  decreased,  and  the  thermal  resistance 
increased  in  the  uninsulated  state  when  compared  to  the  insulated  case.  The  uninsulated 
compensation  chamber  operated  at  a  temperature  3.6°C  lower  than  the  insulated  case  as 
expected  due  to  free  convection  and  radiative  heat  loss.  To  estimate  the  amount  of  heat 
loss  to  the  ambient  in  the  uninsulated  case,  the  exposed  compensation  chamber  was 
modeled  as  a  cylinder  in  free  convection  with  radiation.  The  average  Nusselt  number  for 
free  convection  was  given  by  (Incropera  and  DeWitt,  2002) 

Nud  -  |0.60  +  [1  +  (0.559/Pr)‘'A6]8/27j  ^  ^ 

where 
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with  air  properties  evaluated  at  the  average  temperature  of  the  freestream  and  the  surface. 
The  average  heat  transfer  coefficient  was  given  by 


h  =  -Nud 


(2.13) 


The  total  heat  loss  per  unit  length  from  the  compensation  chamber  was  given  by 

qlot  =  q'conv  +  q'rad  =  hnD(Ts  -  T^o)  +  £nDa(Ts^  -  Tstirr)  (2.14) 
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The  emissivity  of  grade  2  titanium  used  in  this  ealculation  was  £  =  0.3  (Boyer  et  ah, 
1994).  The  heat  loss  from  the  eompensation  chamber  for  the  uninsulated  case  was  found 
to  be  Qcc  =  -6.2  W.  When  the  uninsulated  case  was  included  with  the  previous 
compensation  chamber  heat  input  data,  it  was  found  that  evaporative  heat  transfer 
coefficient  and  thermal  resistance  followed  the  trends  shown  in  Figure  2.14(c).  These 
results  were  expected  since  removing  the  insulation  from  the  compensation  chamber  in 
effect  provided  cooling,  which  moved  the  2(p-lq)  point  away  from  the  condenser.  The 
average  evaporator  and  cold  plate  temperatures  were  significantly  different  which  was 
likely  due  a  10°C  higher  ambient  temperature  during  the  test  involving  the  variation  of 
compensation  chamber  heat  input  from  Qcc  =  0  to  50  W.  As  a  result,  for  this  particular 
case,  it  was  advantageous  to  operate  the  LHP  compensation  chamber  insulated  for 
improved  performance.  For  controlling  the  temperature  of  the  compensation  chamber, 
the  evaporative  heat  transfer  coefficient,  thermal  resistance,  and  operating  temperatures 
were  nearly  identical  between  simultaneous  compensation  chamber  and  evaporator  heat 
input  startup  and  compensation  chamber  temperature  preconditioning,  demonstrating  that 
the  startup  procedure  had  no  impact  on  steady  state  conditions.  However, 
preconditioning  the  compensation  chamber  required  approximately  one  hour  less  time  to 
reach  steady  state  conditions  over  the  simultaneous  heat  input  startup. 

Figure  2.16  presents  transient  LHP  temperatures  for  a  typical  test  at  elevated 
acceleration  (uz  =  1.0  g,  ar  >  0  g).  With  the  LHP  at  ambient  conditions,  the  recirculating 
chiller  in  the  low  temperature  loop  was  set  to  Teg  =  35°C.  Heat  was  applied  as  a  step 
function  (in  this  case,  Qin  =  600  W)  while  simultaneously  starting  the  pump  for  the  high- 
temperature  loop  (w^p  =  0.0077  kg/s).  In  addition,  the  radial  acceleration  was  increased 

to  Ur  =  0.1  g,  which  was  a  nominally  small  value  to  prevent  damage  to  the  power  slip 
rings  (Figure  2.16(a)).  The  LHP  was  allowed  to  achieve  steady  state  conditions  at  = 
0.1  g,  indicated  by  dT/dt  (Figure  2.16(c))  decreasing  to  below  the  threshold  of  0.01 
K/min,  then  the  acceleration  was  increased  to  the  next  desired  radial  acceleration  value 
(in  this  case,  Ur  =  10.0  g).  The  LHP  was  again  allowed  to  achieve  steady  state  conditions 
at  the  given  acceleration  (Figure  2.16(b)),  then  the  acceleration  was  reduced  back  to  ar  = 
0.1  g  for  a  minimum  of  thirty  min.  If  another  elevated  acceleration  was  desired,  steady 
state  at  ar  =  0.1  g  was  reached  before  increasing  the  acceleration  level.  When  the 
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acceleration  was  inereased  to  Ur  =  10.0  g  at  t  =  15,000  s  in  Figure  2.16(b),  the  average 
evaporator  temperature  inereased  by  11°C.  The  2(p-l(p  point  moved  to  the  eondenser 
outlet  from  the  liquid  line  with  increasing  aeceleration,  indicated  by  the  small  oseillations 
in  temperature  at  the  TC13.  The  amount  of  subeooling  inereased  overall  as  indieated  by 
the  deerease  in  temperature  at  the  bayonet  inlet  (TCI 4).  These  phenomena  may  be  due  in 
part  to  fluid  redistribution  in  the  LHP  and  is  discussed  in  the  following  paragraph. 

As  the  rotational  veloeity  of  the  eentrifuge  inereased,  the  resultant  aeeeleration 
vector  magnitude  and  direetion  ehanged  (Figure  2.17(a))  whieh  influeneed  the 
distribution  of  fluid  in  the  LHP.  Subeooled  liquid  entering  the  primary  wiek  of  the 
evaporator  was  forced  to  the  outboard  side  of  the  evaporator  body,  opposite  of  the  heat 
souree,  and  perhaps  leading  to  a  partial  dry-out  of  the  wick  (Figure  2.17(b)).  The 
elevated  aeeeleration  also  hindered  the  ability  of  the  secondary  wick  in  the  eompensation 
ehamber  to  supply  the  evaporator  with  liquid  due  to  pooling.  In  the  eondenser,  pooling 
oeeurred  in  the  bends  of  the  eondenser  eoil,  again  due  to  the  aeeeleration  gradient. 
Depending  on  the  aeeeleration  veetor  direetion,  this  pooling  eould  either  open  or  elose 
the  passage  to  vapor  flow  (Figure  2.17(e)).  All  of  these  phenomena  are  a  result  of 
centrifuge  testing.  Due  to  the  short  radius,  strong  aeeeleration  gradients  oeeur  that  eould 
have  advantageous  or  adverse  effeets  on  the  LHP  operation.  Operation  in  an  aireraft 
environment,  with  signifieantly  larger  radii  during  turns,  will  provide  a  more  uniform 
aeeeleration  gradient  aeross  the  LHP  and  potentially  yield  different  temperature  profiles, 
evaporative  heat  transfer  eoeffieients,  and  thermal  resistanees. 

Figure  2.18  shows  the  thermal  performanee  of  the  LHP  for  radial  aeeelerations 
ranging  from  ar=0.1tol0.0g  and  heat  inputs  ranging  from  Qi^  =  100  to  600  W.  The 
evaporative  heat  transfer  eoeffieient.  Figure  2.18(a),  again  deereased  with  transported 
heat,  similar  to  the  trend  in  Figure  2.13(a)  for  the  stationary  LHP.  The  thermal  resistanee 
of  the  LHP  (Figure  2.18(b)),  was  found  to  deerease  to  a  minimum,  then  inerease,  again 
similar  to  the  stationary  test  results  shown  in  Figure  2.13(b).  In  faet,  when  eombining  the 
stationary  and  elevated  aeeeleration  test  data,  it  was  found  that  the  evaporative  heat 
transfer  eoeffieient  and  thermal  resistance  data  were  in  elose  agreement  with  each  other, 
regardless  of  the  radial  aeeeleration.  This  indicated  that  beneh  top  testing  of  the  LHP  was 
a  reliable  method  for  determining  the  evaporative  heat  transfer  eoeffieient  and  thermal 
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resistance  of  a  LHP  in  an  elevated  acceleration  environment.  However,  it  will  be  shown 
that  this  was  not  true  with  respect  to  finding  the  dry-out  limit.  The  wall  superheat  (Figure 
2.18(c))  was  higher  at  elevated  accelerations  when  compared  to  ar  =  0  g.  This  was 
possibly  due  to  fluid  redistribution  in  the  evaporator  forcing  liquid  away  from  the  heater. 
A  summary  of  all  elevated  acceleration  steady  state  data  points  can  be  seen  in  Table  2.8. 

Ku  et  al.  (2000a  and  2000b)  observed  that  radial  acceleration  changed  the  fluid 
distribution  throughout  the  LHP  which  changed  operating  temperatures  and  that 
acceleration  could  either  increase  or  decrease  LHP  operating  temperatures.  In  addition, 
temperature  overshoots  were  observed  for  mid-range  heat  inputs  and  the  wall  superheat, 
defined  as  the  temperature  difference  between  the  evaporator  and  compensation  chamber, 
was  independent  of  heat  input  and  acceleration  during  startup.  In  the  present  elevated 
acceleration  tests,  it  was  observed  that  the  accelerating  force  changed  the  fluid 
distribution  within  the  LHP,  causing  the  operating  temperatures  to  change.  However,  in 
all  instances,  it  was  observed  that  elevated  acceleration  forces  increased  operating 
temperatures  over  those  at  ar  =  0.1  g.  Significant  temperature  overshoots  were  not 
observed  in  any  of  the  elevated  acceleration  tests. 

Figure  2.19  shows  the  transient  response  of  the  LHP  during  a  series  of  dry-out 
events.  Dry-out  was  indicated  by  a  steady  increase  in  the  evaporator  temperature  and  a 
decrease  in  the  heat  extracted  by  the  calorimeter  Qout-  In  addition,  the  position  of  the  2(p- 
I(p  point  in  the  condenser  moved  toward  the  evaporator  as  indicated  by  a  sequential 
decrease  in  the  condenser  temperatures.  This  occurred  because  the  evaporator  no  longer 
generated  a  sufficient  flow  of  vapor  which  changed  the  operating  point  of  the  LHP.  In 
Figure  2.19(a),  the  LHP  reached  steady  state  while  rotating  slowly  at  ar  =  0.1  g  and  Qm  = 
400  W.  The  rotational  speed  of  the  centrifuge  table  was  increased  until  the  radial 
acceleration  reached  ar  =  8.0  g  at  t  =  300  s.  After  the  evaporator  temperature  TC06 
reached  Ikmax  =  150°C,  the  radial  acceleration  was  reduced  back  to  ar  =  0.1  g.  At  this 
time,  the  evaporator  temperature  continued  to  increase,  but  then  leveled  off  and  then 
decreased  back  to  nearly  the  same  temperature  as  the  previous  steady  state.  In  fact,  all  of 
the  LHP  temperatures  returned  to  within  1°C  of  the  original  steady  state  except  for  TCI  3 
(condenser  outlet),  which  returned  to  within  4°C  of  the  previous  steady  state.  This  larger 
temperature  difference  in  TCI 3  was  attributable  to  a  slight  change  in  the  location  of  the 
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2(p-l(p  point  in  the  condenser.  This  recovery  behavior  shows  that  the  LHP  was  capable  of 
repriming  at  the  end  of  an  acceleration  burst  even  if  the  heat  input  remained  constant.  In 
Figure  2.19(b)  and  Figure  2.19(c),  the  radial  acceleration  was  again  increased  from  Qr  = 
0.1  to  8.0  g  with  the  same  heat  input  (gm  =  400  W).  In  fact,  all  of  the  experiments 
presented  in  Figure  2.19  were  performed  sequentially.  In  Figure  2.19(b),  the  evaporator 
temperature  TC06  was  allowed  to  reach  Te^max  =  175°C  before  decreasing  the  radial 
acceleration  to  ar  =  0.1  g,  and  in  Figure  2.19(c),  the  evaporator  temperature  TC06  reached 
T'e,max  =  200°C  bcforc  decelerating.  In  each  instance,  the  evaporator  temperature 
continued  to  increase,  reached  a  maximum,  and  then  decreased  to  the  original  steady 
state.  However,  the  intensity  of  dry-out  did  seem  to  have  an  impact  on  the  ability  of  the 
LHP  to  reprime.  In  Figure  2.19(c),  with  a  maximum  evaporator  temperature  at 
deceleration  of  Te^max  =  200°C,  the  evaporator  temperature  reached  two  maximums  before 
finally  decreasing  back  to  the  previous  steady  state,  whereas  in  Figure  2.19(a)  and  Figure 
2.19(b),  the  maximum  evaporator  temperatures  reached  a  peak  and  then  monotonically 
decreased.  This  indicated  that  if  the  evaporator  temperature  were  much  higher  than 
200°C,  the  LHP  may  not  have  recovered,  which  would  have  required  that  the  heat  input 
be  reduced  to  zero. 

Figure  2.20  shows  the  temperature  traces  associated  with  the  test  at  =  200  W 
and  ar  =  0.1  and  Qr  =  4.0  g.  Following  the  previously  mentioned  startup  procedures,  the 
LHP  reached  a  quasi-steady  state  while  the  centrifuge  table  rotated  slowly  for  ar  =  0.1  g, 
as  shown  in  Figure  2.20(a).  Similar  to  the  stationary  case  at  this  heat  input,  the  LHP 
temperatures  oscillated,  showing  that  the  heat  pipe  was  operating  during  reversals  in  the 
liquid  flow  due  to  the  liquid-vapor  meniscus  in  the  secondary  wick  moving  back  and 
forth.  Overall,  the  temperatures  shown  in  Figure  2.20(a)  were  quite  close  to  the  case 
shown  in  Figure  2.11(b),  as  presented  in  Table  2.9.  In  addition,  the  period  of  the 
oscillation  of  the  Or  =  0.1  g  case  was  nearly  identical  to  the  ar  =  0  g  case  (approximately 
175  s).  The  only  significant  differences  in  the  independent  variables  between  the  two 
tests  were  the  ambient  temperature  (ATamb  =  5.3°C),  and  the  relatively  small  value  of  the 
radial  acceleration.  Of  note,  however,  was  the  location  of  the  2(p-l(p  point  in  the 
condenser:  For  ar  =  0  g,  this  point  resided  close  to  TC12,  whereas  for  the  case  in  which 
ar  =  0.1  g,  the  2(p-l(p  point  was  near  TCIO.  The  linear  distance  between  these  two  points 
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was  approximately  143  cm.  While  it  was  impossible  to  know  the  exact  location  of  the 
2(p-lq)  point  due  to  the  coarse  resolution  of  the  thermocouples  in  the  condenser,  it  was 
obvious  that  the  location  had  changed  significantly  between  the  two  cases.  In  addition, 
the  evaporative  heat  transfer  coefficient  decreased  and  the  thermal  resistance  increased 
from  ar  =  0  to  0.1  g.  It  was  believed  that  this  small  value  of  the  radial  acceleration 
resulted  in  a  significant  change  in  location  of  the  2(p-l(p  point  due  to  the  pooling  of  liquid. 

After  achieving  the  quasi-steady  state  at  ar  =  0.1  g,  the  radial  acceleration  was 
increased  to  Or  =  4.0  g,  and  the  LHP  again  reached  a  quasi-steady  state,  as  shown  in 
Figure  2.20(b).  The  average  evaporator  temperature  increased  by  more  than  30°C,  and 
the  temperature  in  the  bayonet  inlet  ranged  from  37  <  Tbayonet  inlet  <  70°C,  which  was  a 
much  larger  range  than  that  for  Ur  =  0.1  g.  Oscillations  were  again  seen  at  this 
acceleration  level,  but  the  period  of  the  oscillations  increased  to  approximately  350  s. 
This  may  be  due  to  the  distance  that  the  meniscus  travelled  within  the  evaporator,  which 
resulted  in  wider  swings  in  the  evaporator  temperatures  and  significant  oscillations  of  the 
cold  plate  outlet  temperature,  which  was  nearly  steady  in  the  Or  =  0.1  g  case. 

Figure  2.21  shows  the  steady  state  performance  map  for  the  LHP  relating  radial 
acceleration  and  heat  transported  for  =  2.0  to  10.0  g  and  Qm  =  100  to  600  W.  It  was 
observed  that  dry-out  conditions  occurred  at  varying  radial  accelerations  for  gin  =  100  to 
400  W.  Dry-out  conditions  were  not  observed  through  ar  =  10.0  g  at  gin  =  500  and  600 
W.  Quasi-steady  state  conditions  were  observed  at  gin  =  200  W  and  Ur  =  4.0  g.  This 
demonstrated  that  bench-top  testing  cannot  be  used  to  determine  the  dry-out  limit  with 
respect  to  elevated  acceleration. 

2.5.  Conclusions 

The  effect  of  changes  in  evaporator  heat  input,  compensation  chamber  heat  input, 
and  radial  acceleration  on  a  titanium-water  loop  heat  pipe  were  investigated  for  gin  =100 
to  600  W,  gcc  =  0  to  50  W,  and  ar  =  0.0  to  10.0  g.  A  transient  temperature  rate  of  change 
method  was  developed  to  ensure  steady  state  had  been  achieved.  For  evaporator  heat 
input  gin  =  100  to  600  W,  it  was  observed  that  the  evaporative  heat  transfer  coefficient 
decreased  monotonically,  thermal  resistance  decreased  to  a  minimum,  then  increased 
over  the  same  range,  and  wall  superheat  monotonically  increased.  Flow  reversal  was 
observed  at  gin  =  200  W  due  to  vapor  bubble  generation  in  the  evaporator. 
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When  examining  the  effect  of  compensation  chamber  heat  input  for  gin  =  500  W, 
it  was  found  that  the  average  evaporator  temperatures  dropped  by  15°C  and  evaporative 
heat  transfer  coefficient  improved  by  68%  with  only  a  3%  increase  in  heat  load.  These 
results  differ  from  Ku  (1999)  in  that  an  improvement  was  observed  for  compensation 
chamber  heat  input  up  to  the  point  where  subcooling  was  occurring  in  the  condenser. 
Flow  reversal  was  observed  starting  at  gcc  =  35  W  due  to  the  increased  pressure  in  the 
compensation  chamber  driving  the  liquid/vapor  meniscus  backwards.  Operating  the  LHP 
compensation  chamber  uninsulated  at  gin  =  500  W  was  found  to  degrade  the  LHP 
performance  for  this  particular  case  and  preconditioning  the  compensation  chamber 
temperature  prior  to  evaporator  heat  input  shortened  the  time  to  steady  state. 

When  examining  the  effect  of  radial  acceleration,  it  was  found  that  dry-out 
conditions  occurred  more  readily  at  lower  heat  inputs  (gin  =  100  to  400  W)  than  at  higher 
heat  inputs  (gin  =  500  to  600  W).  The  LHP  was  found  to  be  able  to  reprime  after  an 
acceleration  event  that  caused  dry-out  without  the  heat  input  being  reduced  to  zero.  It 
was  also  observed  that  radial  acceleration  had  little  effect  on  the  evaporative  heat  transfer 
coefficient  and  thermal  resistance  of  the  LHP.  Wall  superheat  was  found  to  be  higher  at 
steady  state  elevated  accelerations  when  compared  to  ar  =  0  g.  This  led  to  conclusion  that 
bench  top  testing  of  the  LHP  is  a  reliable  method  for  determining  the  evaporative  heat 
transfer  coefficient  and  thermal  resistance  of  a  LHP  in  an  elevated  acceleration 
environment  induced  by  a  centrifuge  table,  but  is  not  sufficient  for  determining  the  wall 
superheat  and  dry-out  limit.  These  results  may  or  may  not  actually  occur  in  an  aircraft 
environment  as  centrifuge  operation  can  induce  artifacts  in  the  data  due  to  the  short 
radius  of  operation. 

2.6,  Future  Work 

Experimentation  in  this  thesis  has  been  conducted  using  strict  regimented 
procedures  for  repeatability  and  to  allow  comparisons  across  data  sets.  Typical  operation 
of  LHPs  is  not  along  regimented  schedules  but  in  transient  environments  where  heat 
sources,  heat  sinks,  and  accelerating  forces  are  varying  with  respect  to  time.  As  such, 
experimentation  should  be  conducted  using  transient  profiles  to  more  closely  mimic 
actual  aircraft  environments.  Experimentation  should  also  be  conducted  with  tighter 
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control  over  the  eompensation  ehamber  temperature  to  enhanee  repeatability  and  tailor 
operation  to  speeific  heat  sourees  and  sinks. 
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Figure  2.1.  Loop  heat  pipe  operation.  Adapted  and  reprinted  with  permission  from 
AIAA  (Hoang  and  Ku,  2003). 
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Figure  2.2.  Evaporator  schematic:  (a)  Side  view;  (b)  Cross-sectional  view.  Adapted  and 
reprinted  with  permission  from  AIAA  (Hoang  and  Ku,  2003). 
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(b) 


Figure  2.5.  Thermocouple  locations  on  the  LHP;  (a)  Locations  of  thermocouples  TC04 
through  TCI 5  across  the  LHP;  (b)  Locations  of  TC04  through  TC07  within  the 
evaporator. 
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(C) 


(d) 


Figure  2.6.  Mounting  of  loop  heat  pipe  to  eentrifuge  table,  front  and  top  views:  (a) 
Evaporator  and  eompensation  chamber:  (b)  Transport  lines;  (c)  Condenser  with  cold 
plate;  (d)  Complete  loop  heat  pipe. 
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(a) 


(c) 


Figure  2.7.  High  temperature  fluid  loop:  (a)  Sehematie;  (b)  Reservoir,  pump,  filter, 
flowmeter,  TC03,  and  liquid/liquid  heat  exchanger;  (c)  Cold  plate,  TCOO,  and  TCOl. 
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Figure  2.8.  Use  of  a  cold-start  test  to  determine  when  steady  state  occurred  for  the 
stationary  LHP  (gi„  =  600  W,  Qcc  =  0  W,  ar  =  0  g,  =  0.0077  kg/s,  =  67.7°C,  T^mh 

=  38.1°C):  (a)  Transient  temperature  traces;  (b)  Transient  rate  of  change  of  temperatures; 
(c)  Transient  thermal  resistance  and  evaporative  heat  transfer  coefficient. 
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Figure  2.9.  Transient  startup  of  the  stationary  LHP  (Qin  =  600  W,  Q^c  =  0  W,  ar  =  0  g, 
=  0.0077  kg/s,  T^p  =  67.7°C,  Lamb  =  38.1°C):  (a)  Initial  startup;  (b)  Complete  startup 
until  steady  state. 
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Figure  2.10.  Transient  temperature  profiles  in  the  condenser  and  bayonet  tube  of  the 


stationary  LHP  (Qcc  =  0  W,  ar  =  0  g,  =  0.0077  kg/s,  36.8  < 


<71.6°C,31.7<ra™b 


<  38.1°C):  (a)  =  100  W;  (b)  =  200  W;  (c)  Q,n  =  300  W;  (d)  =  400  W;  (e)  Q,n  = 

500  W;  (f)  gin  =  600  W. 
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Figure  2.1 1.  Transient  temperature  profiles  of  the  stationary  LHP  for  =  200  W  {Qcc  = 
0  W,  Or  =  0  g,  =  0.0077  kg/s,  =  46.1°C,  Lamb  =  31.7°C):  (a)  Transient 

temperature  profiles;  (b)  2(p-l(p  point  oscillation  in  the  condenser. 
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Figure  2.12.  Steady  state  temperature  distribution  versus  transported  heat  for  the 
stationary  LHP  (Qcc  =  0  W,  ar  =  0  g,  =  0.0077  kg/s,  36.8  <  T^p  <  71.6°C,  31.7  <  Lamb 
<  38.1°C):  (a)  Evaporator  seetion;  (b)  Condenser  seetion. 
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Figure  2.13.  Steady  state  performanee  eharaeteristies  of  the  stationary  LHP  versus 
transported  heat  (Qcc  =  0  W,  ar  =  0  g,  =  0.0077  kg/s,  36.8  <  <  67.7°C,  27.6  <  T^mh 

<  38.7°C):  (a)  Evaporative  heat  transfer  coefficient;  (b)  Thermal  resistance;  (c)  Wall 
superheat 
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Figure  2.14.  Steady  state  performance  characteristics  of  the  stationary  LHP  versus 
compensation  chamber  heat  input  (Qin  =  500  W,  =  0  g,  =  0.0077  kg/s,  63.4  <  < 

64.8°C,  36.1  <  Lamb  <  38.1°C):  (a)  Evaporator  temperatures;  (b)  Condenser  temperatures; 
(c)  Evaporative  heat  transfer  coefficient  and  thermal  resistance. 
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Figure  2.15.  Transient  temperature  profiles  in  the  eondenser  and  bayonet  tube  of  the 
stationary  LHP  for  Qcc  =  25  to  50  W  (Qm  =  500  W,  ar  =  0  g,  =  0.0077  kg/s,  63.4  < 

<  64.8°C,  36.1  <  ra„,b  <  38.1°C):  (a)  Q,,  =  25  W;  (b)  Q,,  =  30  W;  (c)  Q,,  =  35  W;  (d) 
0CC  =  40  W;  (e)  Q,,  =  45  W;  (f)  Q,,  =  50  W. 
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Figure  2.16.  Transient  temperature  traces  of  the  LHP  at  elevated  acceleration  (Qi^  =  600 
W,  gee  =  0  W,  =  0.0077  kg/s,  55.2  <  <  59.7°C,  27.9  <  <  30.1°C):  (a)  a,  = 

0.1  g  startup  phase;  (b)  Transition  to  and  steady  state  at  ar  =  10.0  g;  (c)  Transient  rate  of 
change  of  temperatures. 
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Figure  2.17.  Effect  of  resultant  acceleration  vector  direction  on  fluid  distribution  within 
the  LHP:  (a)  Resultant  acceleration  vector  orientation  versus  radial  acceleration;  (b) 
Liquid  pooling  in  the  evaporator,  compensation  chamber,  and  condenser  under  elevated 
acceleration  (to  scale,  top  view);  (c)  Liquid  pooling  in  the  condenser  bends. 
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Figure  2.18.  Steady  state  performanee  characteristics  of  the  LHP  versus  transported  heat 
at  stationary  and  elevated  acceleration  (Q^c  =  0  W,  =  0.0077  kg/s,  37.2  <  < 

67.TC,  25.1  <  Lamb  <  38.7°C):  (a)  Evaporative  heat  transfer  coefficient;  (b)  Thermal 
resistance;  (c)  Wall  superheat. 
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Figure  2.19.  Transient  temperature  traces  of  the  LHP  at  elevated  acceleration  showing 
dry-out  behavior  (Qin  =  400  W,  Qcc  =  0  W,  =  0.0077  kg/s,  312  <  <  59.7°C,  T^mh 
=  28.0°C):  (a)  T„  =  150°C;  (b)  =  175°C;  (c)  T„  =  200°C. 
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Figure  2.20.  Quasi-steady  state  temperature  traees  of  the  LHP  and  eold  plate  at  elevated 
aeeeleration  for  Qin  =  200  W  (Qcc  =  0  W,  =  0.0077  kg/s,  =  41.9°C,  Lamb  = 

26.4°C):  (a)  Transient  temperature  traee  at  Ur  =  0.1  g  and  t  =  13834  s;  (b)  Transient 
temperature  traee  at  =  4.0  g  and  t  =  31240  s. 
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Figure  2.21.  Steady  state  performanee  map  of  the  LHP  relating  radial  aeeeleration  and 
heat  transported  (Q^c  =  0  W,  =  0.0077  kg/s,  37.2  <  <  59.7°C,  25.1  <  Lamb  < 

30.2°C). 
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Table  2.1.  AFRL/RZPS  design  requirements. 


Requirement 

Parameter 

Thermal 

Minimum  Heat  Load 

500  W 

Minimum  Heat  Flux 

rw7cn? 

Operating  Temperature 

200°C 

Condenser  Heat  Sink  Temperature 

5  to  140°C 

Tilt  in  One  G 

±  0  inches,  horizontal 

Conductance 

50°C/W 

Proof  of  Pressure  Test 

3102  psi  (200°C) 

Materials 

Evaporator  Envelope  Material 

Titanium,  CP  Grade  2 

Evaporator  Wick  Material 

Titanium,  CP  Grade  2 

Transport  Eine  Material 

Titanium,  CP  Grade  2 

Working  Eluid 

Water 

EHP  Dimensions 

Evaporator  Configuration 

2.54  cm  OD  up  to  25.4  cm  long 

Evaporator  Eootprint 

20.32  X  10.16  cm 

Condenser  Eootprint 

30.48  X  28.58  cm 

Transport  Eine  Eengths 

Approx.  243.8  cm 
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Table  2.2.  ACT  LHP  geometric  design  parameters. 


Transport  Lines 

Vapor  Line  Length 

Approx.  243.8  cm 

Vapor  Line  Diameter 

0.9525  OD  X  0.0889  cm  wail 

Liquid  Line  Length 

Approx.  335.3  cm 

Liquid  Line  Diameter 

0.6350  OD  X  0.0889  cm  wail 

Condenser  Line  Length 

Approx.  279.4  cm 

Condenser  Line  Diameter 

0.9525  OD  X  0.0889  cm  wall 

Compensation  Chamber 

Diameter 

6.033  cm  OD 

Length 

1 1.43  cm 

Chamber  Location 

Coaxial  with  evaporator 

Wick  Properties 

Effective  Pore  Radius 

9.1  pm 

Permeability 

1.2xl0-'^m^ 

Outside  Diameter 

2.286  cm 

Length 

20.32  cm 

Inside  Diameter 

0.8001  cm 

Grooves 

6 

Groove  Depth 

0.1524  cm 

Groove  Width 

0.1524  cm 
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Table  2.3.  Summary  of  LHP  thermocouple  locations 


Thermocouple 

Location 

TC04 

Evaporator  1 

TC05 

Evaporator  2 

TC06 

Evaporator  3 

TC07 

Evaporator  4 

TC08 

Evaporator  Outlet 

TC09 

Condenser  Inlet 

TCIO 

Condenser  1 

TCll 

Condenser  2 

TC12 

Condenser  3 

TC13 

Condenser  Outlet 

TC14 

Bayonet  Inlet 

TC15 

Compensation  Chamber  /  Evaporator 

65 


Table  2.4.  Summary  of  uncertainties. 


Quantity 

Uncertainty 

Measured 

Temperature 

±  (fixed  value  +  confidence  interval)  °C 

Acceleration 

+  (0.01  +  confidence  interval  +  (ar,oe  -  «r,ci))  g 
-  (0.01  +  confidence  interval  +  (ar,ci  -  «r,ie))  g 

Voltages 

±  (0.00025V  +  0.005  +  confidence  interval)  V 

Constant 

Resistors 

±  0.02%R 

Wick  D,  L 

±  0.0000254  m 

Calculated 

T 

cp 

±  max(Arout,  Afin) 

±  max(A(TC04),  A(TC05),  A(TC06),  A(TC07)) 

ATsh 

±  1.28°C 

Cp,PAO 

±  0.5%  •  c,  (f, )+  (c,  (?;,  +  at;  J  -  c,  (f, )) 

Qoul 

Equation  (2.8) 

h 

Equation  (2.9) 

R 

Equation  (2.10) 
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Table  2.5.  Steady  state  operating  eharaeteristies  for  the  stationary  LHP  (Qcc  =  0  W,  ar  =  0  g,  =  0.0077  kg/s,  36.8  <  <  71.6°C, 

27.6  <  Lamb  ^  38.7°C)  showing  effeet  of  startup  path. _ _ _ _ _ 


Final  gin 
(W) 

Path  to  Final 
gm 

K  (°c) 

T’e/cc  (°C) 

7;p  (°c) 

gout(W) 

k  (W/m^-K) 

7?(K/W) 

100 

0-100 

66.6  ±0.34 

61.8  ±0.86 

39.5  ±0.12 

79.1  ±4.3 

1670  ±390 

0.342  ±0.019 

0-100 

64.9  ±  0.34 

59.9  ±0.86 

37.4  ±0.12 

78.4  ±4.3 

1660  ±380 

0.351  ±0.020 

0-100 

64.8  ±0.34 

59.8  ±0.86 

37.6  ±0.12 

77.9  ±4.3 

1670  ±390 

0.348  ±  0.020 

0-100 

65.4  ±0.34 

60.1  ±0.86 

36.8  ±0.12 

79.0  ±4.2 

1540  ±280 

0.361  ±0.020 

0-100 

65.4  ±0.34 

60.2  ±0.86 

38.2  ±0.12 

77.9  ±4.3 

1540  ±300 

0.350  ±0.020 

0-100-300-100 

66.9  ±0.34 

62.2  ±0.86 

38.8  ±0.12 

79.1  ±4.3 

1730  ±560 

0.355  ±0.020 

0-100-500-100 

67.4  ±  0.34 

62.6  ±0.86 

39.0  ±0.12 

78.8  ±4.3 

1620  ±500 

0.361  ±0.020 

133 

0-100-133 

65.1  ±0.34 

59.0  ±0.86 

40.4  ±0.12 

114±5.5 

1770  ±310 

0.217  ±0.011 

166 

0-100-133-166 

65.6  ±0.34 

58.3  ±0.86 

43.0  ±0.12 

143  ±6.5 

1760  ±250 

0.158  ±0.0076 

200 

0-200 

66.3  ±  0.46 

66.3  ±0.86 

46.1  ±0.12 

174  ±7.7 

1560  ±230 

0.116  ±0.0058 

0-200 

66.4  ±0.35 

66.4  ±0.86 

45.7  ±0.12 

174  ±7.7 

1460 ± 140 

0.119  ±0.0057 

0-200 

65.5  ±0.34 

65.5  ±0.86 

45.4  ±0.12 

176  ±7.8 

1660  ±250 

0.114  ±0.0055 
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Table  2.5,  continued.  Steady  state  operating  characteristics  for  the  stationary  LHP  (Qcc  =  0  W,  ar  =  0  g,  =  0.0077  kg/s,  36.8  < 

<  71.6°C,  27.6  <  Tamh  <  38.7°C)  showing  effect  of  startup  path.  _ _ _ _ 


Final  gi„ 
(W) 

Path  to  Final 

Qm 

T,  (°C) 

T’e/cc  (°C) 

T.,  (°C) 

0out(W) 

h  (W/m^-K) 

7?(K/W) 

300 

0-100-300 

77.1  ±0.39 

62.4  ±0.86 

51.6±0.12 

260  ±  1 1 

1660 ± 160 

0.0985  ±  0.0045 

0-100-300 

75.8  ±0.40 

60.9  ±0.86 

50.2  ±0.12 

259  ±  11 

1670 ± 170 

0.0988  ±  0.0045 

0-200-300 

78.1  ±0.38 

63.1  ±0.86 

52.0  ±0.12 

258  ±  11 

1580 ± 140 

0.101  ±0.0046 

0-300 

79.3  ±  0.40 

64.4  ±0.86 

52.8  ±0.12 

258  ±  11 

1550 ± 160 

0.103  ±0.0047 

400 

0-200-400 

94.0  ±  0.34 

73.6  ±0.86 

57.8  ±0.12 

345  ±  14 

1400  ±  85 

0.105  ±0.0045 

0-300-400 

94.3  ±0.35 

75.0  ±0.86 

58.4  ±0.12 

345  ±  14 

1470  ±  98 

0.104  ±0.0045 

0-400 

97.4  ±  0.34 

72.6  ±0.86 

56.5  ±0.12 

344  ±  14 

1190  ±60 

0.119  ±0.0051 

500 

0-100-300-500 

116±0.34 

82.4  ±0.86 

61.7  ±0.12 

430  ±  18 

1050  ±45 

0.126  ±0.0053 

0-100-500 

119±0.34 

83.8  ±0.86 

63.1  ±0.12 

429  ±  18 

999  ±41 

0.131  ±0.0055 

0-400-500 

117±0.34 

83.6  ±0.86 

62.4  ±0.12 

431  ±  18 

1050  ±47 

0.127  ±0.0053 

0-500 

122  ±0.34 

85.3  ±0.86 

63.7  ±0.12 

432  ±  18 

956  ±  42 

0.135  ±0.0057 

600 

0-200-400-600 

145  ±0.36 

93.1  ±0.86 

67.8  ±0.12 

515±21 

803  ±  33 

0.149  ±0.0062 

0-400-500-600 

141  ±0.34 

93.3  ±0.86 

67.5  ±0.12 

515±21 

867  ±  37 

0.142  ±0.0060 

0-600 

148  ±0.36 

93.7  ±0.86 

67.7  ±0.12 

513±21 

757  ±33 

0.157  ±0.0066 
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Table  2.6.  Steady  state  operating  characteristics  for  the  stationary  LHP  showing  effect  of  heat  input  to  the  compensation  chamber  (Qi^ 


=  500  W,  ar=  0  g,  =  0.0077  kg/s,  63.4  <  <  64.8°C,  36.1  <  <  38.1°C). 


Qoo  (W) 

t;  (°c) 

Te/cc  (°C) 

7;p  (°c) 

Qont  (W) 

k  (WW-K) 

R(K/W) 

0 

122  ±0.34 

85.3  ±0.86 

63.7  ±0.12 

432  ±  18 

956  ±  42 

0.135  ±0.0057 

5 

115  ±0.35 

85.3  ±0.86 

63.6  ±0.12 

438  ±  18 

1130±58 

0.118  ±0.0050 

10 

109  ±0.35 

85.6  ±0.86 

63.8  ±0.12 

444  ±  18 

1420  ±  no 

0.103  ±0.0044 

15 

107  ±0.35 

85.1  ±0.86 

63.4  ±0.12 

449  ±  19 

1610±  130 

0.0970  ±0.0041 

20 

108  ±0.34 

86.1  ±0.86 

64.0  ±0.12 

456  ±  19 

1610±  130 

0.0972  ±0.0041 

25 

no  ±0.34 

87.7  ±0.86 

64.2  ±0.12 

461  ±  19 

1620 ± 130 

0.0991  ±0.0042 

30 

112  ±0.34 

90.0  ±0.86 

64.5  ±0.12 

466  ±  19 

1610±  130 

0.103  ±0.0044 

35 

118±0.35 

95.5  ±0.86 

64.6  ±0.12 

468  ±  19 

1580 ± 120 

0.114  ±0.0048 

40 

120  ±0.34 

97.4  ±0.86 

64.8  ±0.12 

473  ±  20 

1580 ± 120 

0.117  ±0.0049 

45 

127  ±0.34 

104  ±0.86 

64.5  ±0.12 

475  ±  20 

i530±  no 

0.132  ±0.0055 

50 

129  ±0.35 

105  ±0.86 

64.5  ±0.12 

477  ±  20 

i5io±  no 

0.134  ±0.0056 
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Table  2.7.  The  effect  of  compensation  chamber  temperature  control  on  LHP  operation  (Qin  =  500  W,  ar=  0  g,  =  0.0077  kg/s,  T^p 
=  52.5°C,  =  26.4°C)  _ 


Compensation  Chamber 
Conditions 

Qco  (W) 

K  CC) 

Te,  CC) 

Te.  CC) 

0out(W) 

k  (W/m^-K) 

T(K/W) 

Time  to 
Steady 
State 

Uninsulated,  no 
temperature  control 

0 

115 

52.4 

59.3 

442  ±  19 

878  ±37 

0.142  ±0.0060 

300  min. 

Insulated,  no 
temperature  control 

0 

107 

52.6 

62.9 

447  ±  19 

1050  ±51 

0.122  ±0.0052 

375  min. 

Insulated,  temperature 
controlled  to  Tec  = 

72.8°C,  simultaneous 
heat  input  startup 

20 

103 

53.5 

72.8 

470  ±  20 

1340  ±77 

0.106  ±0.0045 

310  min. 

Insulated,  temperature 
controlled  to  Tec  = 

72.8°C,  preconditioned 
CC 

100 

decreased 
to  20 

103 

53.4 

72.8 

467  ±  19 

1350  ±77 

0.106  ±0.0045 

250  min. 
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Table  2.8.  Steady  state  operating  eharacteristies  of  the  rotating  LHP  {Qcc  =  0  W,  =  0.0077  kg/s,  37.2  <  <  59.7°C,  25.1  <  Lamb 

<30.2°C).  _ _ _ _ _ _ _ 


an(W) 

«r  (g) 

K  (°c) 

Te/cc  (°C) 

7;p  (°c) 

gout  (W) 

k  (W/m^-K) 

7?(K/W) 

100 

0.1  ±0.012 

64.9  ±  0.34 

59.8  ±0.86 

38.0  ±0.12 

74.2  ±  4.2 

1620  ±340 

0.362  ±0.021 

100 

0.1  ±0.012 

64.6  ±  0.34 

59.8  ±0.86 

38.1  ±0.12 

74.6  ±4.2 

1810±510 

0.355  ±0.020 

100 

2.0  ±0.050 

66.4  ±  0.34 

61.2  ±0.86 

37.5  ±0.12 

79.9  ±4.5 

1770  ±420 

0.361  ±0.021 

100 

4.1  ±0.091 

71.2  ±0.44 

64.3  ±0.86 

37.2  ±0.16 

76.1  ±4.8 

1130±  150 

0.446  ±  0.029 

200 

0.1  ±0.012 

70.8  ±0.41 

59.3  ±0.86 

42.0  ±0.13 

174  ±7.8 

1400  ±  no 

0.165  ±0.0078 

200 

0.1  ±0.012 

71.1  ±0.43 

59.5  ±0.86 

41.9±0.13 

172  ±7.8 

i380±  no 

0.169  ±0.0081 

200 

0.1  ±0.012 

71.5  ±0.38 

60.3  ±0.86 

42.0  ±0.12 

174  ±7.8 

1360±  no 

0.170  ±0.0079 

200 

2.0  ±0.050 

75.8  ±0.65 

61.6  ±0.86 

41.7  ±0.21 

188  ±9.2 

1220  ±86 

0.181  ±0.0096 

300 

0.1  ±0.012 

76.4  ±  0.48 

57.2  ±0.86 

46.0  ±0.12 

275  ±  12 

1390 ± 100 

0.111  ±0.0051 

300 

2.0  ±0.050 

81.7  ±0.76 

61.2  ±0.86 

45.6  ±0.16 

285  ±  13 

1280  ±80 

0.126  ±0.0063 

400 

0.1  ±0.012 

86.4  ±0.39 

60.6  ±0.86 

49.3  ±0.12 

360  ±  15 

1300  ±82 

0.103  ±0.0045 

400 

0.1  ±0.012 

86.6  ±0.39 

61.3  ±0.86 

49.3  ±0.12 

361  ±  15 

1300  ±81 

0.103  ±0.0045 

400 

0.1  ±0.012 

86.1  ±0.39 

61.1  ±0.86 

49.3  ±0.12 

361  ±  15 

1330  ±84 

0.102  ±0.0044 

400 

0.1  ±0.012 

85.8  ±0.42 

61.1  ±0.86 

49.3  ±0.12 

361  ±  15 

1330  ±89 

0.101  ±0.0044 

400 

2.0  ±0.050 

87.6  ±  0.42 

59.7  ±0.86 

49.3  ±0.12 

372  ±  16 

1280  ±74 

0.103  ±0.0045 

400 

4.0  ±0.089 

91.5  ±0.43 

59.8  ±0.86 

49.3  ±0.12 

376  ±  16 

1160  ±56 

0.112  ±0.0051 

400 

4.1  ±0.091 

93.0  ±0.63 

60.0  ±0.86 

49.3  ±0.15 

376  ±  17 

n30±59 

0.116  ±0.0055 

400 

6.0  ±0.13 

94.8  ±0.34 

60.7  ±0.86 

49.9  ±0.12 

369  ±  16 

1030  ±46 

0.122  ±0.0053 
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Table  2.8,  continued.  Steady  state  operating  characteristics  of  the  rotating  LHP  {Qcc  =  0  W,  =  0.0077  kg/s,  37.2  <  <  59.7°C, 

25.1  <ra^b<30.2°C). _ 


an(W) 

«r  (g) 

T,  (°C) 

Te/cc  (°C) 

T.,  (°C) 

Qout  (W) 

h  (W/m^-K) 

7?(K/W) 

500 

0.1  ±0.012 

110±0.36 

69.3  ±0.86 

52.9  ±0.12 

440  ±  18 

960  ±  43 

0.131  ±0.0055 

500 

0.1  ±0.012 

107  ±0.35 

69.0  ±0.86 

52.6  ±0.12 

446  ±  19 

1030  ±48 

0.123  ±0.0052 

500 

0.1  ±0.012 

107  ±0.41 

69.2  ±0.86 

52.7  ±0.12 

447  ±  19 

1030  ±50 

0.122  ±0.0052 

500 

0.1  ±0.012 

109  ±0.42 

69.5  ±0.86 

52.9  ±0.12 

447  ±  19 

1000  ±46 

0.125  ±0.0053 

500 

0.1  ±0.012 

108  ±0.35 

68.9  ±0.86 

52.6  ±0.12 

448  ±  19 

1010  ±47 

0.123  ±0.0052 

500 

2.0  ±0.050 

114  ±0.42 

65.5  ±0.86 

52.5  ±0.12 

452  ±  19 

887  ±37 

0.136  ±0.0058 

500 

4.0  ±  0.090 

110±0.34 

64.4  ±  0.86 

53.1  ±0.13 

459  ±  19 

982  ±  43 

0.123  ±0.0053 

500 

6.0  ±0.13 

111  ±0.34 

64.6  ±0.86 

54.0  ±0.13 

460  ±  20 

934  ±  42 

0.125  ±0.0055 

500 

8.0  ±0.17 

115  ±0.34 

65.1  ±0.86 

54.8  ±0.12 

463  ±  20 

865  ±38 

0.131  ±0.0058 

500 

10  ±0.21 

117±0.34 

64.7  ±0.86 

55.2  ±0.12 

452  ±  20 

824  ±  35 

0.136  ±0.0060 

600 

0.1  ±0.012 

137  ±0.36 

76.6  ±0.86 

56.0  ±0.12 

533  ±22 

759  ±32 

0.152  ±0.0063 

600 

0.1  ±0.012 

141  ±0.43 

76.6  ±0.86 

55.9  ±0.12 

533  ±22 

707  ±  29 

0.160  ±0.0067 

600 

0.1  ±0.012 

141  ±0.43 

76.3  ±0.86 

55.8  ±0.12 

533  ±22 

703  ±  29 

0.160  ±0.0067 

600 

0.1  ±0.012 

140  ±0.42 

76.2  ±0.86 

55.8  ±0.12 

531  ±22 

707  ±  30 

0.159  ±0.0067 

600 

0.1  ±0.012 

133  ±0.34 

76.2  ±0.86 

56.0  ±0.12 

539  ±23 

821  ±35 

0.143  ±0.0060 

600 

2.0  ±0.050 

140  ±0.34 

72.9  ±0.86 

56.3  ±0.12 

547  ±  24 

726  ±31 

0.152  ±0.0066 

600 

4.0  ±  0.090 

152  ±0.50 

71.7  ±0.86 

57.1  ±0.12 

546  ±  23 

621  ±26 

0.173  ±0.0075 

600 

6.0  ±0.13 

148  ±0.34 

71.5  ±0.86 

57.8  ±0.12 

545  ±  24 

651  ±27 

0.165  ±0.0072 

600 

8.1  ±0.17 

147  ±0.34 

71.8  ±0.86 

58.8  ±0.12 

546  ±  24 

662  ±  28 

0.162  ±0.0072 

600 

10  ±0.21 

145  ±0.34 

72.2  ±  0.86 

59.7  ±0.13 

538  ±24 

691  ±29 

0.158  ±0.0070 
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APPENDIX  A,  OPERATING  PROCEDURES 


A.l,  Standard  Operating  Procedure 

0.  Ensure  the  main  power  breaker  is  in  the  “OFF”  position.  The  breaker  is  located 
on  the  71B  H-bay  second  floor  landing.  See  Figure  A.l  for  picture.  Ensure  the  control 
panel  box  is  in  “Man”  mode,  potentiometer  is  set  to  zero  (turned  completely  counter¬ 
clockwise),  and  the  E-stop  button  has  been  depressed.  The  control  panel  box  is  located  in 
Control  Room  2.  See  Figure  A.2  for  picture  of  control  panel  box. 

1 .  Mount  test  payload  with  sufficient  factor  of  safety  to  ensure  a  reliable  mounting 
configuration  consistent  with  the  generated  forces.  The  retaining  method  will  be 
designed  for  a  payload  mass  subjected  to  a  maximum  table  capacity  of  12  g. 

2.  Mount  the  appropriate  counter  balance  weight  at  the  appropriate  location  to  place 
the  centroid  in  the  center  of  the  table  and  level  each  spoke  to  minimize  vibrational  noise 
and  table  runout. 

3.  Ensure  all  centrifuge  maintenance  has  been  completed. 

4.  Software  Startup  Procedure 

a.  Open  the  Fab  VIEW  VI  file  needed  to  control  the  data  acquisition  and 
table  voltage. 

b.  It  is  up  to  the  experimentalist  to  write  the  program  used  to  control  data 
acquisition  and  table  voltage.  Sample  programs  that  perform  these  tasks 
are  available.  See  Figure  A.3  for  sample. 

c.  Press  the  ‘SYSTEM  ENABEED’  button  on  the  front  panel  so  that  it  is 
illuminated. 

d.  Press  the  ‘Run’  button  located  in  the  top-left  of  the  EabVIEW  toolbar. 
This  will  cause  the  program  to  become  functional. 

5.  Prior  to  each  set  of  experimental  testing: 

a.  Check  all  bolts  to  ensure  all  experimental  apparatus  are  tightened 
properly. 
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b.  Clear  the  table  of  foreign  objeets,  ineluding  removing  all  tools  and  plaeing 
in  the  appropriate  CTK.  Ensure  all  free  standing  equipment  and  furniture 
are  seeurely  plaeed  along  the  perimeter  of  the  room. 

c.  Turn  table  power  switeh  on.  This  provides  power  to  the  instrumentation 
and  deviees  on  the  table.  See  Figure  A.4  for  loeation. 

d.  Cheek  out  instrumentation  for  proper  operation. 

e.  Turn  on  reeireulating  ehiller  loeated  on  the  east  wall.  See  Figure  A.5  for 
pieture.  Cheek  ehill  bath  eoolant  flows  and  flow  rates  of  any  intermediate 
flow  loops.  A  booster  pump  for  the  ehill  bath  eoolant  is  available.  Follow 
the  chill  bath  plumbing  schematic  in  Figure  A.  6  to  enable  operation.  The 
booster  pump  control  is  remotely  located  in  Control  Room  2.  See  Figure 
A. 7  for  picture. 

f  Fock  inter-connect  door. 

g.  Check  camera  operation. 

h.  Make  final  check  on  the  table  for  tools  or  loose  objects. 

6.  Unlock  main  power  breaker  and  flip  to  ‘ON’  position.  See  Figure  A.l  for 
location. 

a.  Ensure  the  emergency  stop  button  is  activated  and  the  potentiometer  is 
turned  completely  counter-clockwise  on  the  control  panel  box  prior  to 
proceeding  with  powering  the  centrifuge  table  motor. 

b.  Turn  on  warning  beacon  and  evacuate  personnel.  Warning  beacon 
switches  are  located  outside  Test  Cell  4  and  on  the  west  wall  of  control 
room  1. 

c.  Engage  table  motor  control  switch  on  the  north  wall  (cooling  motor  will 
be  operational).  See  Figure  A.8  for  location. 

d.  Secure  outer  doors. 

e.  Place  “Test  in  Progress  Do  Not  Enter”  sign  on  the  outer  door. 

7.  Controlling  the  System 

a.  To  control  the  voltage  supplied  to  the  table,  turn  the  mode  switch  from 
‘MAN’  to  ‘AUTO’  on  the  control  panel  box  and  flip  the  ‘Table  Voltage’ 
switch  to  the  ‘ON’  position.  While  the  ‘Table  Voltage’  switch  is  ‘ON’, 
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the  voltage  ean  be  adjusted  using  the  slider  bar  on  the  left  side  called 
‘Table  Voltage’.  Table  voltage  can  be  turned  on  and  off  as  often  as 
desired  while  the  system  is  engaged. 

8.  Conduct  test  procedure.  Test  procedures  are  experiment  dependent  and  up  to  the 
experimentalist  to  develop. 

9.  Shutdown 

a.  Slowly  reduce  the  table  voltage  to  zero  using  the  slider  bar.  Flip  the 
‘Table  Voltage’  switch  to  the  ‘OFF’  position  on  the  Lab  VIEW  VI. 

b.  On  the  control  panel  box,  turn  the  switch  from  ‘AUTO’  to  ‘MAN’  and 
press  the  ‘STOP’  button.  Wait  for  the  table  to  stop  rotating  then  press  the 
‘E-STOP’  button. 

c.  Press  the  ‘SYSTEM  ENABLED’  button  on  the  EabVIEW  VI  so  that  it  is 
no  longer  illuminated.  The  program  will  stop. 

10.  Emergency  Shutdown 

a.  If  for  any  reason  an  emergency  should  occur  press  the  ‘E-STOP’  on  the 
control  panel  box.  Should  the  table  “run  away”  or  suddenly  accelerate  the 
motor  will  automatically  shutdown.  Contact  the  appropriate  personnel 
prior  to  a  restart  after  an  emergency  shutdown. 

•  Andrew  Eleming  58942 

•  Earry  Byrd  53238 

•  Travis  Michalak  64429 

A,2.  Test  Procedures 

1 .  Stationary  Operation 

a.  Set  chiller  to  T^g  =  35°C.  Allow  to  come  to  steady  state. 

b.  Turn  high  temperature  coolant  loop  on  and  set  to  m  =  0.0077  kg/s. 
Simultaneously,  apply  desired  heat  load  to  evaporator. 

c.  Allow  EHP  to  achieve  steady  state  operation  by  examining  dT/dt  plot  of 
TCOO,  TCOl,  TC04,  TC05,  TC06,  TC07,  TC08,  and  TC09.  Steady  state  is 
achieved  when  all  of  these  thermocouples  are  bracketed  by  -0.01  <  dT/dt  < 
0.01  K/min. 
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d.  Remove  heat  load  from  the  evaporator  and  turn  off  the  high  temperature 
eoolant  loop  or  adjust  to  next  desired  heat  load  and  repeat  (e). 

2.  Rotational  Operation 

a.  Set  ehiller  to  T^g  =  35°C.  Allow  to  eome  to  steady  state. 

b.  Inerease  radial  aeeeleration  to  Or  =  0. 1  g. 

e.  Turn  high  temperature  eoolant  loop  on  and  set  to  m  =  0.0077  kg/s. 
Simultaneously,  apply  desired  heat  load  to  evaporator. 

d.  Allow  LHP  to  aehieve  steady  state  operation  by  examining  dT/dt  plot  of 
TCOO,  TCOl,  TC04,  TC05,  TC06,  TC07,  TC08,  and  TC09.  Steady  state  is 
aehieved  when  all  of  these  thermoeouples  are  braeketed  by  -0.01  <  dT/dt  < 
0.01  K/min. 

e.  Inerease  radial  aeeeleration  to  desired  level. 

f  Again,  allow  LHP  to  aehieve  steady  state  operation  by  examining  dT/dt 
plot  of  TCOO,  TCOl,  TC04,  TC05,  TC06,  TC07,  TC08,  and  TC09.  Steady 
state  is  aehieved  when  all  of  these  thermoeouples  are  braeketed  by  -0.01  < 
dT/dt  <  0.01  K/min. 

g.  Deerease  radial  aeeeleration  to  ar  =  0. 1  g. 

h.  If  shutting  down,  allow  LHP  to  operate  for  30  min.,  then  remove  heat  load 
from  the  evaporator  and  turn  off  the  high  temperature  eoolant.  Deerease 
radial  aeeeleration  to  Ur  =  0.0  g.  If  eontinuing  testing,  allow  LHP  to 
aehieve  steady  state  operation  by  examining  dT/dt  plot  of  TCOO,  TCOl, 
TC04,  TC05,  TC06,  TC07,  TC08,  and  TC09.  Steady  state  is  aehieved 
when  all  of  these  thermoeouples  are  braeketed  by  -0.01  <  dT/dt  <  0.01 
K/min.  Repeat  (e)-(h). 
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(b) 

Figure  A.l.  Centrifuge  table  main  power  breaker:  (a)  Eleetrical  panel  MCC-6;  (b) 
Centrifuge  table  main  power  breaker. 
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Figure  A.2.  Centrifuge  table  control  panel  box. 
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Figure  A.3.  Sample  Lab  VIEW  eontrol  program. 
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Figure  A.4.  Centrifuge  table  power  switeh. 
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Figure  A.5.  Neslab  recirculating  chiller. 
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To  Centrifuge  Table 


From  Centrifuge  Table 


Figure  A.6.  Chill  bath  plumbing  schematic. 
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Figure  A.7.  Booster  pump  eontrol  panel. 
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Figure  A.8.  Centrifuge  table  motor  control  power  switch. 


88 


APPENDIX  B,  UNCERTAINTY  ANALYSIS 


The  uncertainty  analysis  for  this  experiment  was  broken  into  the  following  three 
different  categories:  measured,  constant,  and  calculated  quantities.  For  each  steady  state 
condition,  151  data  points  were  collected  representing  five  minutes  of  data.  Measured 
values  were  averaged  and  uncertainties  were  calculated  based  on  the  fixed  error  of  each 
instrument  plus  the  confidence  interval  for  the  average.  The  fixed  error  for 
thermocouples  and  mass  flow  meter  was  determined  via  the  calibration  methods  and  can 
be  seen  in  Appendix  C.  The  fixed  error  of  the  accelerometer  was  ±0.1  g  and  ±0.00025*V 
+  0.005  volts  for  voltage  measurements.  The  confidence  interval  was  based  off  a 
statistical  t-distribution  with  confidence  level  of  0.95  and  was  given  by  Montgomery 
(Montgomery,  2003) 


Confidence  Interval  =  ±t  (B.l) 

where  t  is  a  tabulated  value  based  on  the  confidence  level  and  number  of  degrees  of 
freedom,  a  is  the  sample  standard  deviation,  and  n  is  the  number  of  data  points  in  the 
sample.  Constant  quantities  included  precision  resistors  used  for  current  measurement 
for  heat  input  via  evaporator,  compensation  chamber,  and  preheater  with  an  uncertainty 
of  ±0.02%  and  wick  diameter  and  length  measurements,  with  an  uncertainty  of 
±0.0000254  m.  Calculated  quantities  used  uncertainty  methods  that  would  be  the  most 
conservative  for  the  experiment.  For  averaged  quantities  including  T^,  and  h  ,  the 

largest  uncertainty  of  the  individual  measurements  was  used  as  the  uncertainty  of  the 
average  value.  For  the  specific  heat  of  PAO,  the  uncertainty  was  0.5%  of  the  total  value 
plus  difference  between  the  specific  heat  using  the  upper  limit  of  the  temperature 
measurement  and  the  average  specific  heat.  The  calorimetry  of  the  cold  plate,  gout,  was 
given  by 


l2out 


^l^p,PAo(^out  ^in) 


(B.2) 
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The  uncertainty  of  gout  was  determined  to  be 


^Qout  — 
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out 
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where 
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The  average  evaporative  heat  transfer  coefficient,  h  ,  was  given  by 

aout 


h  = 


nDL(_T,  -  T^) 


The  uncertainty  of  h  was  given  by 


(B.8) 
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(B.ll) 


_  Qout 

^  “  TiD^Lif^  -  T^) 

dh  _  -Qout 
dL  nDL^(T^-T^) 

_  Qout 

W,  ~  tzdl(t,  -  T^y 

dh  _  -Qout 

dT^  nDL(f,  -  T^y 

The  thermal  resistanee,  R,  was  given  by 

f  —f 

^  ^cp 

l2out 

The  uneertainty  of  the  thermal  resistanee  was  given  by 


where 

dR  _-(Te-fcp) 
5  Qout  QLt 

dR  _  1 

dR  _  -1 
dfcp  Qout 


(B.12) 
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APPENDIX  C.  CALIBRATION  OF  THERMOCOUPLES  AND  FLOW  METER 


C.l,  Thermocouple  Calibration 

The  calibrations  of  the  thermocouples  used  on  the  centrifuge  table  required  new 
programming  since  the  data  acquisition  system  was  upgraded  for  the  centrifuge  table. 
Control  of  a  calibration  bath  and  RTD  were  integrated  with  the  centrifuge  table’s  data 
acquisition  unit.  A  Lab  VIEW  program  was  generated  to  simultaneously  control  all  three 
devices. 

The  thermocouples  were  calibrated  using  a  Hart  Scientific  6330  Calibration  Bath 
and  Hart  Scientific  1502A  NIST-Traceable  platinum  resistance  temperature  detector 
(RTD).  The  bath  was  capable  of  producing  steady  state  temperatures  from  40  to  280°C. 
The  calibration  bath  used  Dow  Corning  200.50  silicon  oil.  There  were  several  steps 
required  before  a  thorough  calibration  of  the  thermocouples  could  be  determined.  First, 
Lab  VIEW  software  needed  to  be  written  to  interface  with  the  calibration  bath,  RTD,  and 
data  acquisition  system.  Second,  the  characteristics  of  the  calibration  bath  and  RTD 
needed  to  be  determined  with  respect  to  response  times  and  temperature  fluctuations. 
Finally,  complete  calibration  curves  for  each  of  the  thermocouples  needed  to  be 
developed. 

In  the  loop  heat  pipe  experimental  setup,  there  were  four  grounded  probe 
thermocouples  and  twelve  exposed  tip  type  E  thermocouples.  The  grounded  probe 
thermocouples  were  used  in  the  calorimetry  of  the  cold  plate,  temperature  measurement 
for  the  flow  meter  calibration,  and  the  ambient  temperature  inside  the  box,  where  the 
error  needed  to  be  reduced  as  much  as  possible.  As  a  result,  it  was  decided  to  only 
calibrate  these  four  over  the  anticipated  temperature  range,  from  20  to  145°C,  in  5°C 
intervals.  The  twelve  exposed  tip  thermocouples  were  mounted  on  the  loop  heat  pipe  in 
various  locations.  They  needed  to  be  calibrated  over  the  full  20  to  230°C  temperature 
range.  Thus,  the  calibration  of  the  thermocouples  needed  a  second  calibration  bath  that 
was  capable  of  achieving  the  temperature  range  of  20  to  40°C.  The  Hart  Scientific  6330 
Calibration  Bath  was  not  capable  of  maintaining  a  steady  temperature  below  40°C  for 
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Dow  Corning  200.50  bath  fluid.  The  bath  used  for  this  portion  of  the  ealibration  was  a 
Brinkmann  Lauda  RCS  20-D  ealibration  bath  using  Brayeo  Mieronie  889 
(polyalphaolefin  or  PAO)  as  the  bath  fluid.  The  same  RTD  was  used  as  before,  and  the 
temperature  set  point  for  the  ealibration  bath  was  set  manually  rather  than  through  the 
data  aequisition  interfaee.  After  some  examination,  it  was  notieed  that  due  to  the  limited 
eapability  of  the  Brinkmann  bath  to  heat  and  eool,  it  was  diffieult  to  maintain  a  fairly 
eonstant  temperature  in  the  PAO.  As  a  result,  an  insulated  eopper  tube  with  a  elosed  base 
was  plaeed  in  the  bath,  and  the  eopper  tube  was  filled  with  PAO.  This  signifieantly 
redueed  the  temperature  fiuetuation  in  the  bath  temperature. 

The  first  step  in  the  development  of  the  thermoeouple  ealibrations  was  to  write 
the  Lab  VIEW  software  to  interfaee  with  the  ealibration  bath,  RTD,  and  data  aequisition. 
Initially,  sub- Vi’s  were  developed  to  interfaee  direetly  with  the  ealibration  bath  for 
reading  the  eurrent  temperature  and  setting  the  bath  temperature  set  point  as  well  as  for 
reading  the  RTD.  Wire  diagrams  of  the  RTD  read,  ealibration  bath  temperature  set  point, 
and  ealibration  bath  eurrent  temperature  read  ean  be  found  in  Figure  C.l  through  Figure 
C.3,  respeetively.  Onee  these  sub-VI’s  were  developed,  they  needed  to  be  ineorporated 
into  a  larger  framework.  The  front  panel  and  wire  diagram  for  the  automatie  ealibration 
program  ean  be  seen  in  Figure  C.4  while  the  front  panel  and  wire  diagram  for  the  manual 
ealibration  program  ean  be  seen  in  Figure  C.5. 

Due  to  eleetrieal  noise  in  the  eentrifuge  table  test  eell  and  fiuetuations  in  the 
ealibration  bath  temperature,  a  method  was  devised  to  determine  that  the  ealibration  bath 
had  reaehed  a  steady  state.  The  previous  100  RTD  temperature  values  were  reeorded  into 
an  array,  and  the  standard  deviation  of  the  sampling  was  ealeulated.  This  number  of 
RTD  readings  was  seleeted  for  a  95%  probability  and  a  eonfidenee  interval  of  0.95 
(Montgomery,  2003).  When  the  standard  deviation  dropped  below  the  speeified 
threshold  of  0.005°C,  the  program  would  indieate  that  the  system  had  reaehed  steady 
state  and  the  thermoeouples  were  then  read.  The  bath  temperature  would  then  be 
ineremented  or  deeremented  as  neeessary.  First,  the  full  ealibration  proeess  began  by 
plaeing  the  thermoeouples  in  the  Hart  Seientifie  ealibration  bath  for  the  temperature 
range  40  to  230°C  in  5°C  inerements  ramping  up  and  down  with  a  standard  deviation 
threshold  of  0.005°C.  Then,  the  thermoeouples  and  RTD  were  eleaned  and  plaeed  in  the 
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Brinkmann  bath  for  the  temperature  range  20  to  35°C.  The  data  from  eaeh  proeess  was 
eombined  to  produee  one  eomposite  data  set  eonstituting  the  entire  temperature  range. 

For  data  reduction,  all  100  data  points  constituting  one  nominal  temperature  value 
were  averaged  and  the  confidence  interval  was  calculated  for  the  RTD.  After  this 
process,  due  to  the  increment  and  decrement  of  the  calibration  process,  there  were  two 
data  points  for  each  nominal  temperature  value,  as  shown  in  Figure  C.6.  Plots  of  RTD 
versus  each  thermocouple  were  generated,  and  polynomial  trend  lines  were  fit  for  each 
thermocouple.  A  sample  plot  of  TCOO  is  given  in  Figure  C.7.  Also,  in  an  effort  to  reduce 
maximum  deviation  of  the  actual  versus  calculated  RTD  values,  higher-order  polynomial 
trend  lines,  from  first  to  fifth  order,  were  implemented  to  evaluate  maximum  deviation. 
A  fifth-order  polynomial  was  selected  since  it  reduced  the  maximum  deviation  by 
approximately  a  factor  of  4  over  a  first-order  trend  line.  Results  from  this  analysis  can  be 
seen  in  Table  C.l  for  TCOO.  Table  C.2  shows  the  trend  line  equations  for  each 
thermocouple  in  tabular  form. 

The  uncertainty  associated  with  each  thermocouple  was  determined  by  accounting 
for  four  sources  of  error:  maximum  measured  uncertainty  inherent  to  the  RTD,  the 
maximum  confidence  interval  of  the  RTD  temperature  over  the  100  readings  used  in  the 
calibration,  the  maximum  deviation  of  the  calculated  temperature  from  the  measured 
temperature,  and  the  confidence  interval  associated  with  the  100  data  points  in  the  sample 
of  the  thermocouples.  The  error  inherent  to  the  RTD  was  ±0.009°C.  The  maximum 
confidence  interval  of  the  RTD  from  the  100  readings  was  ±0.0055°C  over  the  entire 
temperature  range.  The  maximum  deviations  of  the  calculated  temperature  and  the 
measured  temperature  are  thermocouple  specific,  and  can  be  seen  in  Table  C.3  along  with 
the  total  error  of  each  thermocouple. 

C.l.l.  Calibration  Procedure 

1.  Mount  all  of  the  thermocouples  to  the  RTD  probe  with  the  thermocouple  and 
RTD  tips  as  close  to  each  other  as  possible. 

2.  Place  the  thermocouple  and  RTD  bundle  vertically  into  the  Hart  Scientific  6330 
Calibration  Bath,  with  the  probes  not  touching  any  of  the  bath  surfaces. 

3.  Turn  on  the  RTD  and  calibration  bath,  setting  the  bath  to  40°C. 

4.  Turn  on  the  centrifuge  table  power. 
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5.  In  the  control  room,  turn  on  Lab  VIEW  and  open  the  thermocouple  calibration 
program. 

6.  Set  the  lower  temperature  set  point  to  40°C  and  the  upper  temperature  set  point  to 
230°C,  with  a  temperature  increment  of  5°C. 

7.  Set  the  standard  deviation  threshold  to  0.005°C. 

8.  Run  the  VI. 

9.  Current  temperatures  can  be  examined  while  the  calibration  is  in  progress  on  the 
main  screen. 

10.  After  this  calibration  cycle  has  been  completed,  remove  and  clean  the 
thermocouples  and  RID,  rebundle,  and  place  in  the  copper  tube  that  is  located  in 
the  Brinkmann  Lauda  RCS  20-D  calibration  bath. 

1 1 .  Turn  the  bath  on  and  set  to  20°C. 

12.  Open  the  manual  thermocouple  calibration  VI. 

13.  Ensure  the  “Keep  Running”  button  is  depressed,  then  run  the  program. 

14.  Wait  for  the  RED  standard  deviation  to  reach  below  0.005,  the  press  the  ‘Proceed 
to  TC  Read’  button. 

15.  Increment  the  bath  temperature  by  5°C  up  to  35°C  and  back  down  to  20°C, 
recording  temperature  values  at  each  location  an  overall  two  times. 

16.  Take  the  two  data  files  and  combine  to  make  one  composite  data  file. 

C.2,  Flow  Meter  Calibration 

The  calibration  of  the  turbine  flow  meter  for  the  high-temperature  fluid  loop 
proved  to  be  a  difficult  challenge.  Due  to  the  chemistry  of  PAO,  there  is  a  significant 
difference  in  density  and  viscosity  with  respect  to  temperature.  As  a  result,  a  calibration 
surface  that  was  dependent  on  flow  meter  output  voltage  and  temperature  was  required  to 
determine  the  actual  mass  flow  rate.  An  uncertainty  analysis  was  also  performed  to  fully 
characterize  the  flow  meter. 

A  calibration  setup  was  developed  using  a  Lauda  RCS-20D  calibration  bath  filled 
with  PAO,  Tuthill  pump  from  the  high-temperature  fluid  loop,  SS-56S6  Whitey  inline 
filter  with  a  140  micron  filter,  calibrated  type-E  thermocouple  that  was  used  for 
measuring  temperature  entering  the  flow  meter  in  the  high  temperature  fluid  loop,  and  an 
Omega  ETB-9506  turbine  flow  meter  with  ELSC-61  signal  conditioner.  A  Lab  VIEW 
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code  was  developed  to  aid  in  the  calibration  of  the  flow  meter  (Figure  C.8).  This  eode 
served  as  the  stopwateh  and  data  reeording  for  the  ealibration.  Temperature  and  flow 
meter  voltage  were  reeorded,  with  raw  data  fdes  generated.  The  pump  was  controlled  via 
a  0-1 OV  input  signal  to  a  TECO  FM50  motor  eontroller.  The  flow  meter  required  an 
input  voltage  between  12-28VDC  and  output  voltage  between  O-IOV.  All  data  was 
eolleeted  through  the  eentrifuge  table  slip  rings  to  the  data  aequisition  to  eapture  all  error 
inherent  to  the  eentrifuge  table  test  bed.  A  ealibrated  Mettler  PC4400  seale  was  used  to 
determine  the  mass  eolleeted  during  a  given  test  run.  Flow  straightening  seetions 
upstream  and  downstream  were  placed  aecording  to  the  manufaeturer's  instruetions.  A 
three-way  valve  was  installed  after  the  flow  meter,  whieh  allowed  the  entire  flow  system 
to  reach  a  steady  temperature.  A  schematie  of  the  setup  is  shown  in  Figure  C.9.  Onee 
the  temperature  was  steady,  the  flow  was  diverted  to  a  catch  basin  for  a  specified  amount 
of  time.  The  voltage  from  the  flow  meter  and  the  temperature  from  the  thermoeouple 
were  reeorded  during  this  time,  and  when  the  basin  was  full,  the  flow  was  again  diverted 
to  reeireulating  the  PAO  baek  to  the  ehiller.  During  eaeh  measurement,  as  many  data 
points  as  possible  were  eolleeted  aeross  the  time  span  with  the  limiting  faetor  being  the 
iteration  time  on  the  Fab  VIEW  software.  The  minimum  number  of  data  points  eolleeted 
for  any  given  run  was  437.  The  voltages  and  temperatures  were  averaged  and  a 
eonfidenee  interval  of  0.95  was  ealculated  based  on  a  statistieal  t-distribution  for  eaeh  test 
run.  The  test  was  repeated  for  a  total  of  five  averaged  data  points  for  each  nominal 
temperature  and  flow  rate.  These  tests  were  completed  over  the  range  of  r=  20  to  120°C 
in  intervals  of  25°C  and  flow  rates  ranging  from  m  =  0.0064  to  0.025  kg/s  in  intervals  of 
approximately  0.002  kg/s.  A  3-D  paraboloid  regression  equation  was  generated  using 
SigmaPlot  to  relate  temperature,  flow  meter  voltage,  and  mass  flow  rate  (Figure  C.IO). 
The  equation  ean  be  seen  in  Table  C.4. 

During  the  eourse  of  developing  the  ealibration  setup  and  testing,  several  issues 
were  encountered.  Immediately  from  the  start  of  working  out  the  bugs  in  the  system,  it 
was  notieed  that  the  flow  meter  would  not  output  a  voltage  linearly  as  expected,  but 
rather  responded  in  a  quadratie  fashion  with  a  loeal  maximum  at  approximately  six  volts. 
It  was  discovered  that  the  motor  housing  of  the  pump  was  not  properly  grounded,  eausing 
eleetromagnetic  interferenee  to  disrupt  the  operation  of  the  flow  meter.  The  output 
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voltage  was  extremely  ineonsistent,  ehanging  +/-  0.5  volts  at  apparently  random  times. 
Initially  an  SS-56S6  Whitey  inline  filter  with  a  140  mieron  filter  was  plaeed  in  the  fiow 
loop  to  eateh  any  debris  that  may  have  been  pieked  up  by  the  pump.  It  was  then  notieed 
that  on  the  fiow  loop  outlet,  air  bubbles  were  being  rejeeted.  All  of  the  fittings  were 
retightened  and  the  flow  remained  steady.  At  apparently  random  times  the  fiow  meter 
would  start  outputting  random  voltage  readings.  It  was  diseovered  that  serew  terminal 
was  not  tightened  down  on  the  wire  eonneetor,  oeeasionally  ereating  an  open  eireuit  loop. 
During  the  eourse  of  data  eolleetion,  the  output  voltage  would  develop  a  trend,  as  seen  in 
Figure  C.l  1.  When  this  was  observed  during  data  reduetion,  another  data  point  would  be 
eolleeted  to  replaee  it.  This  was  likely  due  to  fiuetuations  in  the  flow  rate  and  was 
observed  more  prevalently  at  higher  flow  rates. 

With  this  type  of  ealibration,  it  is  eritieal  to  have  a  firm  grasp  on  the  uneertainty 
assoeiated  with  the  mass  fiow  rate.  Three  types  of  uneertainty  were  identified  assoeiated 
with  this  ealibration:  error  assoeiated  with  the  seale  and  time,  error  due  to  the  voltage 
eonfidenee  interval  and  thermoeouple  error  and  eonfidenee  interval,  and  deviation  of  the 
fit  equation  from  aetual  data.  Eaeh  of  these  errors  are  deseribed  below. 

The  total  error  assoeiated  with  this  uneertainty  analysis  is  given  by 

Arhtot  =  Arh^/t  +  Arhy/T  +  Arhdev  (C.l) 

The  error  for  Am^^  was  determined  by 


(C.2) 


(C.3) 


and 


^^m/t  _  1 

dm  t 


(C.4) 
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(C.5) 


drh^/t  _  _rn 
dt 


The  error  for  Arhy/j  was  determined  by  performing  an  error  analysis  on  the  fit  equation 
developed  using  SigmaPlot 


rh  =  Yo  +  aT  +  bV  +  cT^  +  dV^ 
whieh  yields  the  uneertainty 


Arhy/T 


+ 


/ drhy/j 
V  dV 


(C.6) 


(C.7) 


where 


and 


dlfly/j 

dT 


=  a  +  2cT 


(C.8) 


dlfly/j 

dV 


=  b  +  2dV 


(C.9) 


The  error  assoeiated  with  the  deviation  of  the  fit  equation  from  the  aetual  data  is  given  by 

Arhdev  =  “  ^p|  (C.IO) 

Applying  this  uncertainty  analysis  to  the  data  set,  a  maximum  error  of  4.0%  was 
calculated  and  imposed  over  the  entire  calibration  range. 


C.2,1,  Calibration  Procedure 

1 .  Connect  monitor,  keyboard,  and  mouse  to  extensions  in  Test  Cell  4. 

2.  Activate  Lab  VIEW  code  entitled  “Flow  meter  Calibration.” 

3.  Ensure  three-way  valve  is  in  by-pass  mode  (flow  diverting  back  to  tank). 

4.  Turn  on  Eauda  calibration  bath  and  set  to  20°C.  Turn  on  the  motor  control  unit. 

5.  Using  the  EabVIEW  program,  increase  the  motor  control  input  voltage  to  the 
desired  setting  (1.0). 

6.  Allow  the  flow  meter  to  come  to  temperature  using  PAO  from  the  calibration  bath 
that  is  by-passed  back  to  the  bath. 
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7.  Take  the  empty  2  L  flask  and  plaee  on  the  seale.  Tare  the  scale.  Remove  the 
flask  from  the  scale  and  place  under  the  three-way  valve. 

8.  When  ready,  quickly  move  the  three-way  valve  into  measurement  mode  while 
simultaneously  depressing  the  “Timer”  button  on  the  Lab  VIEW  front  panel. 

9.  Allow  the  flask  to  fill  for  45  seconds  for  pump  voltage  settings  one  to  five,  and 
for  30  seconds  for  pump  voltages  settings  six  to  ten. 

10.  Quickly  return  the  three-way  valve  to  by-pass  mode  while  simultaneously 
depressing  the  “Timer”  button  on  the  Lab  VIEW  front  panel. 

1 1 .  User  will  be  prompted  to  enter  the  mass  collected.  Carefully  place  the  filled  flask 
on  the  scale  and  record  the  mass  in  the  program.  Return  the  collected  PAO  to  the 
calibration  bath. 

12.  Repeat  steps  6  through  12,  increasing  the  motor  control  input  voltage  by 
increments  of  l.OV  up  to  10. OV  for  a  given  temperature,  then  increasing  the 
temperature  by  25°C  up  to  120°C. 
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IConFigure  Serial  port  (baud  rate^  I 

Idata  bits,  parity ,  stop  bite  and  Bow  conljoi) .  I 

Figure  C.l.  Lab  VIEW  sub-VI  wire  diagram  for  RED  read. 


Enable  Termination  Char  (T) 


Configure  Serial  port  (baud  rate, 

data  bits,  parity,  stop  bits  and  flow  control). 


Figure  C.2.  Lab  VIEW  sub-VI  wire  diagram  for  calibration  bath  temperature  set. 
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Enable  Termination  Char  (T) 

EH 


Figure  C.3.  Lab  VIEW  sub-VI  wire  diagram  for  calibration  bath  temperature  read. 
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(b) 

Figure  C.4.  LabVIEW  VI  for  controlling  the  automatic  thermocouple  calibration:  (a) 
Front  panel;  (b)  Wire  diagram 
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(a) 


(b) 

Figure  C.5.  LabVIEW  VI  for  manual  thermocouple  calibration;  (a)  Front  panel;  (b) 
Wire  diagram. 
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Figure  C.6.  RTD  temperature  vs.  time  from  the  thermoeouple  calibration  procedure. 
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Figure  C.7.  Sample  RID  vs.  thermocouple  plot  for  TCOO. 
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Figure  C.8.  Lab  VIEW  VI  for  flow  meter  ealibration  program;  (a)  Front  panel;  (b)  Wire 
diagram. 
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From  DAQ 


To  DAQ 


Figure  C.9.  Schematic  of  flow  meter  calibration  loop 
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I 

•B 


T(C) 


V(V) 


Figure  C.IO.  Temperature  and  flow  meter  voltage  versus  mass  flow  rate  ealibration 
curve  for  the  high-temperature  fluid  loop  flow  meter. 
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Figure  C.ll.  Sample  data  eolleeted  during  one  time  run  for  the  flow  meter  ealibration. 
(a)  “Shotgun  Blast”  good  data  set;  (b)  “Trend”  bad  data  set. 
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Table  C.l.  Maximum  deviation  of  ealculated  RTD  and  experimental  RTD  eorresponding 
to  eaeh  order  of  polynomial  for  thermoeouple  TCOO. 


Polynomial  Order 

First 

Seeond 

Third 

Fourth 

Fifth 

Maximum  Deviation  (°C) 

0.48 

0.45 

0.16 

0.15 

0.11 
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Table  C.2.  Coefficients  for  the  trend  line  of  each  thermocouple. 
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Table  C.3.  Maximum  deviation  and  total  error  for  each  thermocouple. 


Thermocouple 

Maximum 

Deviation 

±Total  Error 

TCOO  (CP  In) 

0.I064I4 

0.122388 

TCOl  (CP  Out) 

0.098267 

0.IIII8I 

TC02  (Box  Ambient) 

0.097236 

0.III085 

TC03  (Flow  Meter  In) 

0.098736 

0.III083 

TC04  (Evap  I) 

0.295217 

0.307238 

TC05  (Evap  2) 

0.280751 

0.293427 

TC06  (Evap  3) 

0.280739 

0.292962 

TC07  (Evap  4) 

0.322196 

0.334702 

TC08  (Evap  Out) 

0.315632 

0.328080 

TC09  (Cond  In) 

0.274506 

0.288721 

TCI0(Cond  I) 

0.278507 

0.291299 

TCI  I  (Cond  2) 

0.278103 

0.290722 

TCI 2  (Cond  3) 

0.256994 

0.269226 

TCI 3  (Cond  Out) 

0.291867 

0.304493 

TCI 4  (Bayonet  In) 

0.285927 

0.298477 

TCI 5  (CC/Evap) 

0.285186 

0.297719 
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Table  C.4.  3-D  paraboloid  regression  equation  for  high-temperature  fluid  loop  flow 
meter. 


m  =  +  aT  +  bV  -\-cT^  +  dV^ 


yo 

2.07704738E-03 

a 

-4.69012732E-05 

b 

2.35000226E-03 

e 

1.916501 17E-07 

d 

1.4981 1559E-05 
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APPENDIX  D.  LOOP  HEAT  PIPE  MOUNTING 

The  loop  heat  pipe  was  mounted  onto  the  eentrifuge  table  sueh  that  the  eenterline 
of  the  tubing  eoineided  with  the  outer  table  radius  as  mueh  as  possible.  A  small  deviation 
existed  sinee  the  eondenser  seetion  and  the  evaporator/eompensation  ehamber  were  both 
straight.  This  indueed  a  non-uniform  radial  aoeeleration  field  over  the  lengths  of  these 
seetions  that  needed  to  be  quantified.  As  shown  in  Figure  D.l,  the  straight  eondenser 
section,  with  length  L^,  was  geometrically  aligned  on  the  table  so  that  the  centerline  was 
as  close  to  the  table  radius,  i?ct,  as  much  as  possible  over  its  length.  The  radius  of  the 
condenser  midpoint,  was  found  as  shown  in  Figure  D.l  as  well.  The  condenser 
endpoints  were  first  set  to  coincide  with  the  centrifuge  table  radius.  The  half-angle  is 
given  by 


The  length  from  the  center  of  the  centrifuge  table  to  the  midpoint  of  the  condenser  is 

Lem  =  Ret  cos  e  (D.2) 

One-half  of  the  change  in  radius  from  this  point  to  the  centrifuge  table  radius  is  used  to 
determine  the  radius  of  the  midway  point  of  the  condenser  section. 

i?cm=ii+^(^ct-icm)  (D-3) 

The  evaporator  section  and  the  compensation  chamber  were  also  straight  and  the  method 
to  locate  these  components  in  relation  to  the  outer  table  radius  is  similar  to  that  described 
above  for  the  condenser  section.  However,  the  evaporator  section  is  directly  attached  to 
the  compensation  chamber  and  the  length  of  the  evaporator  is  different  than  that  of  the 
compensation  chamber.  Therefore,  further  care  was  taken  in  determining  the  maximum 
deviation  of  the  centerline  radius  of  these  two  components  from  the  radius  of  the 
centrifuge  table. 
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Stands  were  designed  using  G-7  phenolie  to  mount  the  loop  heat  pipe  vertieally 
with  support  at  the  compensation  chamber,  evaporator,  condenser,  and  two  on  the 
transport  lines  (Figure  D.2).  The  tops  of  these  stands  were  anchored  to  the  table  to 
reduce  deflection  when  the  table  was  rotating.  A  survey  was  taken  at  22  locations  on  the 
loop  heat  pipe  to  determine  how  far  various  portions  of  the  loop  heat  pipe  were  from  the 
centerline  radius  and  at  what  distance  that  location  was  from  the  origin  at  point  1  (Figure 
2.6(d)).  The  loop  heat  pipe  had  a  minimum  radius  of  1 19.2  cm  at  locations  5  and  6  and 
maximum  radius  of  123.3  cm  at  location  15  at  the  outside  edge  of  the  tubing  at  each 
location.  The  entire  loop  heat  pipe  fitted  within  4.6  cm  for  a  percent  radial  difference  of 
3 .7%.  Complete  survey  data  can  be  seen  in  Table  D.  1 . 
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Figure  D.2.  LHP  survey  locations. 
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APPENDIX  E,  BRAYCO  MICRONIC  889  TECHNICAL  DATA 


Technical  data  for  various  properties  against  temperature  for  Brayco  Micronic 
889,  including  density,  specific  heat  at  constant  pressure,  thermal  conductivity,  and 
kinematic  viscosity  were  provided  by  Ghajar  et  al  (1994).  Brayco  Micronic  889  technical 
data  from  Castrol  was  compiled  and  the  data  was  curve-fitted  using  a  least  squares 
approach  and  varying  order  polynomials.  The  equations  are  only  valid  on  the 
temperature  range  -54  <  T  <  135^C.  The  authors  report  that  for  p,  Cp,  and  k  a  maximum 
deviation  of  less  than  ±0.5%  from  the  reported  data  for  the  property  equations.  For  v, 
they  reported  a  maximum  deviation  of  ±3.4%.  The  equations  for  these  properties  are 
given  as 

p  =  1.36  •  10^  -  4.56r  ±  0.01577^  -  0.280  •  lO-^T^ 


±0.174-10-^7^ 

(E.l) 

k  =  0.154-  5.88-  10-^7 

(E.2) 

Cp  =  1.022  ±  3.77  -  10-^7 

(E.3) 

For  all  of  these  equations,  T  is  in  Kelvin,  p  is  in  kg/m^,  Cp  is  in  kJ/kg-K,  and  k  is  in  W/m- 
K. 
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APPENDIX  F,  CENTRIFUGE  TABLE  UPGRADES 


The  previous  centrifuge  table  data  acquisition  system  dated  back  to  the  early 
1990’s.  The  original  computer  was  a  Pentium  386  running  ViewDAC  for  data  collection 
and  reduction.  At  the  start  of  this  project,  it  was  determined  that  the  data  acquisition 
system  needed  to  be  upgraded,  including  the  data  acquisition  unit  and  computer  system. 
Since  all  of  this  equipment  was  going  to  be  upgraded,  it  was  decided  to  completely 
evaluate  the  existing  data  acquisition  wiring  and  document  this  information. 

The  first  step  to  updating  the  data  acquisition  system  was  to  record  the  original 
data  acquisition  wiring.  All  of  the  wiring  on  the  centrifuge  table  was  rewired  and 
documented  so  that  it  would  be  easier  to  trace  wiring  back  to  the  centrifuge  table  control 
room.  The  wiring  scheme  on  the  centrifuge  table  can  be  seen  in  Figure  F.  1 .  The  terminal 
strip  at  the  bottom  of  this  photo  is  located  on  the  rotating  table,  while  the  terminal  strips 
at  the  top  are  on  the  stationary  support.  These  are  connected  via  a  40-ring  slip  ring.  A 
new  wiring  panel  was  developed  for  the  wiring  coming  from  the  stationary  terminal  strips 
above  the  centrifuge  table  back  to  the  control  room,  as  shown  in  Figure  F.2.  Wiring  on 
the  centrifuge  table  terminal  strip  now  matches  the  wiring  coming  into  the  centrifuge 
table  control  room.  This  information  was  completely  documented  for  future  reference. 
Each  circuit  from  the  centrifuge  table  to  the  slip  ring  wiring  panel  in  the  control  room 
was  checked  for  continuity  by  hand. 

After  verifying  the  wiring  configuration  was  in  proper  operating  condition, 
documentation  for  the  new  data  acquisition  hardware  was  reviewed  so  that  the  upgrades 
could  be  started.  Initially,  a  new  wiring  panel  was  developed  for  the  new  data  acquisition 
interface  as  shown  in  Figure  F.3.  This  panel  was  designed  to  accommodate  64  channels 
(three  wires  per  channel)  for  data  acquisition,  as  well  as  16  channels  (four  wires  per 
channel)  for  voltage  and  current  control.  Each  circuit  from  the  data  acquisition  interface 
panel  to  the  two  new  data  acquisition  modules  was  checked  for  continuity  by  hand. 
Jumper  cables  were  created  to  transfer  signals  coming  in  from  the  slip  rings  to  the  data 
acquisition  wiring  panel. 
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The  new  data  aequisition  system  from  Agilent  Teehnologies  has  a  mainframe 
(E8408A)  with  four  slots  into  which  the  following  cards  were  installed;  a  command 
module  (E1406A),  an  8/16  channel  D/A  converter  (E1418A),  a  SVi-digit  multimeter 
(E1411B),  and  a  sixty-four  channel,  3-wire  multiplexer  (E1476A).  The  command 
module  serves  as  the  main  source  of  communication  between  all  of  the  cards  in  the 
mainframe.  The  command  module  also  exchanges  data  and  commands  between  the 
computer  and  data  acquisition  system.  The  D/A  converter  is  a  control  type  module, 
allowing  the  user  to  request  a  certain  output  voltage  or  current  from  up  to  sixteen 
channels.  The  multimeter  reads  the  voltages  that  the  multiplexer  collects  from 
thermocouples,  pressure  transducers,  accelerometers,  etc.,  as  well  as  any  externally 
applied  voltages,  currents,  and  resistances.  Communication  between  the  data  acquisition 
system  and  the  computer  takes  place  via  the  general  purpose  interface  bus  (GPIB)  or 
IEEE-488  protocol.  Essentially,  text  commands  are  sent  from  the  computer  to  the  data 
acquisition  system.  Then,  if  the  command  requests  a  control  signal,  the  proper  output  is 
processed.  If  the  command  were  for  data  acquisition,  then  ASCII  data  is  returned  to  the 
user  for  processing. 

With  the  new  data  acquisition  system  and  computer  assembled,  documentation  for 
the  data  acquisition  system  needed  to  be  reviewed  to  determine  the  proper  commands 
necessary  to  use  the  computer  to  communicate  with  the  data  acquisition  system.  Initially, 
single  text  lines  were  sent  from  the  computer  using  Agilent’s  VISA  Assistant  software. 
Commands  were  sent  to  read  thermocouple  temperatures  on  one  of  the  channels  of  the 
multiplexer. 

The  software  used  for  writing  the  data  acquisition  code  was  EabVIEW,  a  visual 
computer  language.  Virtual  instruments  (Vis)  for  communicating  with  the  data 
acquisition  system  started  fairly  crudely.  Virtual  instruments  are  subprograms  that  are 
written  with  certain  inputs  and  calculated  outputs,  which  can  greatly  simplify  a 
complicated  code.  Eirst,  the  task  of  reading  several  voltage  channels  and  outputting  the 
data  to  the  screen  was  accomplished.  Second,  the  reading  of  several  voltage  channels 
was  placed  inside  a  timed  loop  such  that  data  would  be  recorded  at  regular  intervals  and 
written  to  a  file  that  Microsoft  Excel  could  read.  Next,  the  proper  conversions  for  voltage 
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to  temperature,  acceleration,  and  pressure  were  written  into  the  code  so  that  actual  data 
was  recorded  to  a  fde. 

The  next  task  to  be  accomplished  with  this  software  was  control  communication. 
What  was  desired  was  a  system  where  a  certain  voltage  could  be  applied  to  increase  the 
angular  velocity  of  the  centrifuge  table,  and  to  control  a  variety  of  other  devices,  such  as 
pumps,  heaters,  etc.  Following  a  similar  process  to  the  development  of  the  data 
acquisition  software,  a  virtual  instrument  was  written  that  could  control  the  output 
voltage  of  several  channels  along  with  the  capability  of  turning  them  on  and  off  at  any 
time.  This  VI  was  then  merged  with  the  data  acquisition  program  with  appropriate 
Boolean  commands  for  control.  This  VI  was  tested  and  verified  when  data  was  collected 
for  a  liquid-vapor  separator  experiment  on  the  centrifuge  table.  A  program  was  tailored 
for  this  experiment,  including  the  appropriate  flow  meter,  pressure  transducer, 
accelerometer  conversions,  and  data  recording. 

After  control  of  the  centrifuge  table  was  accomplished  using  a  voltage  from  the 
D/A  converter  to  control  the  angular  velocity,  it  was  decided  to  control  the  acceleration 
directly.  A  relation  was  then  developed  between  voltage  and  acceleration.  Centrifugal 
acceleration  is  given  by 


2 

ro) 

Ur  =  - 

a 

Voltage  is  related  to  angular  velocity  by 


V  ^  Bo) 


(F.l) 


(F.2) 


where  B  is  an  experimental  constant  to  be  determined.  Substituting  this  relation  into  the 
expression  for  centrifugal  acceleration  yields 


(F.3) 


Solving  equation  (F.3)  for  voltage  yields  the  relation  used  for  deriving  a  corresponding 
voltage  for  a  chosen  centrifugal  acceleration. 


(F.4) 
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Thus,  voltage  is  linearly  related  to  the  square  root  of  eentrifugal  aeeeleration.  This 
relation  ean  be  extended  to  relate  voltage  and  the  aeeeleration  magnitude.  The  magnitude 
of  acceleration  is  given  by 


a 


+ 


(F.5) 


where  ,  and  ajare  the  accelerations  normalized  by  gravity  in  the  azimuthal, 

radial,  and  axial  directions  on  the  centrifuge  table.  Solving  for  the  radial  acceleration 
gives 

(F.6) 


When  the  centrifuge  table  is  rotating  with  a  constant  velocity,  =  0  and  =  -1.  After 

substituting  these  values  in  for  equation  (F.6)  and  solving  for  the  square  root  of 
centrifugal  acceleration 

-  1)'/^  (F.7) 

Substituting  equation  (F.7)  into  (F.4)  yields  the  relation  between  table  voltage  and  the 
magnitude  of  the  acceleration. 


To  calculate  B,  experimental  data  relating  table  voltage  with  centrifugal  acceleration  was 
collected.  Then,  a  plot  of  voltage  versus  the  square  root  of  centrifugal  acceleration  was 
generated  and  a  linear  best  fit  regression  was  derived  with  the  voltage  intercept  forced 

through  the  origin.  The  corresponding  slope  is  then  B^Jg  /  r  .  A  sample  plot  can  be  seen 

in  Figure  F.4.  After  the  slope  for  this  equation  was  found,  an  acceleration  control  slide 
bar  was  added  to  the  data  acquisition  program  by  deriving  a  corresponding  voltage 
output.  It  is  important  to  note  that  the  slope  is  experiment  specific,  and  if  the  location  of 
the  accelerometer  is  changed,  a  new  slope  needs  to  be  found. 

It  was  desired  to  have  the  capability  of  reading  higher  temperatures  on  the 
centrifuge  table.  The  current  thermocouple  amplifier  on  the  centrifuge  table  is  for  Type 
T  thermocouples,  which  have  an  operating  temperature  range  between  -250  to  350°C.  A 
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new  Type  E  thermoeouple  amplifier  (Omega  OM7-47-E-07-2-C)  has  been  installed  on 
the  underside  of  the  centrifuge  table  opposite  of  the  existing  thermocouple  amplifier  so 
that  either  one  can  be  used,  depending  on  the  experimental  requirements  specified, 
providing  operating  temperatures  between  -200  to  900°C. 
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I 


Figure  F.l.  Updated  wiring  on  the  eentrifuge  table. 
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Figure  F.2.  Wiring  panel  from  centrifuge  table  to  the  centrifuge  table  control  room. 
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Figure  F.3.  Wiring  panel  for  the  new  data  aequisition  system. 
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3.5 
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APPENDIX  G.  LABVIEW  PROGRAMS 


130 


(b) 

Figure  G.l.  Lab  VIEW  VI  for  the  LHP  experiment;  (a)  Front  panel;  (b)  Wire  diagram. 


I3I 


(a) 


(b) 

Figure  G.2.  Lab  VIEW  sub-VI  wire  diagram  for  voltage  output  eontrol;  (a)  Output  on; 
(b)  Output  off 
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Figure  G.3.  Lab  VIEW  sub-VI  wire  diagram  for  data  aequisition  eommunication. 
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Figure  G.4.  Lab  VIEW  sub-VI  wire  diagram  for  data  analyzing. 


^^educed  Datal 
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Figure  G.5.  Lab  VIEW  sub-VI  wire  diagram  for  data  reeording. 
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APPENDIX  H.  CENTRIFUGE  WIRING  TABLES 
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Table  H.2.  Data  acqusition  terminal  board  wiring. 
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APPENDIX  I,  SAMPLE  CALCULATIONS 


Example  Calculation  of  Average  Heat  Transfer  Coefficient  for  Flat  Plate  Flow 
Given; 

Altitude; 

Mach  number; 

Wall  temperature 
Plate  length; 

Calculations; 

Freestream  temperature; 

Ta,=  (7.7664E-4)//-  (2.8994E-2)/7^  +  (5.3483E-1)//^  -  (9.5033)77+  (4.85  07E+1) 
7’«,=  284.37  K 

Freestream  density; 

/>oo  =  (-4.9336E-6)77^+  (2.0898E-3)77^  -  (8.9917E-2)77+  1.0868 
yOoo  =  0.6870  kg/m^ 

Film  temperature; 

f  =  FooCO.S  +  0.039  Maoo^)  +  (284.37  K)(0.5+0.039(0.  8)V0.5(408.15  K) 

f  =  353.36  K 

Air  density  at  the  film  temperature; 

p*=pJ(Tjf)  =  (0.6870  kgW)(284.37  K)/(353.36  K)  =  0.5528  kgW 
Freestream  speed  of  sound; 

a^=^{yRT^)  =  V((1.4)(286.9  mW-K){  284.37  K))  =  337.97  m/s 
Freestream  velocity; 

7/00  =  Macaco  =  (0.8)(337.97  m/s)  =  270.37  m/s 

Freestream  absolute  viscosity  (Reference  values  from  Incropera  and  DeWitt,  2002); 
pa.  =  pRiTJTKf^^  =  (184.6  X  10'^  N-s/m^)(284.37  K  /  300  K)°-’'’ 
p^=  1.772  X  lO'^N-s/m^ 


77=  5  km 

Maoo  =  0.8 

7k=135°C  =  408.15  K 

7.  =  1.0  m 


Reynolds  number; 

ReL=  (poo7/ooF)//roo=  (0.6870  kg/m^)(270.37  m/s)(l  m)  /  (1.772  x  10'^  N-  s/m^) 
ReL=  1.05  X  10’ 

[TURBUFENT] 
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Prandtl  number  at  the  film  temperature; 

Pr*  =  (-1.2593E-9)7^  +  (1.7778E-6)7^  -  (9.4177E-4)r+  (8.6418E-1) 

Pr*  =  0.6978 

Specific  heat  at  the  film  temperature; 

Cp*  =  (4.4444E-7)7^  -  (3.3333E-5)7^  -  (6.9921E-2)r+  (1.0187E+3) 

Cp*  =  1009.44  J/(kg-K) 

Recovery  factor  for  turbulent  flow; 
r  =  Pr*'^^  =  (0.6978)^^^  =  0.8870 

Adiabatic  wall  temperature; 

Eaw  =  roo[l+r((7-l)/2)Maoo^]  =  (284.37  K)(l+(0.8870)(0.4/2)(0.8f )  =  316.66  K 
Absolute  viscosity  at  the  film  temperature; 

H*  =  =  (184.6  X  lO'^N-  sW)(353.36  K  /  300  =  2.091  x  10'^  N-sW 

Eocal  skin  friction  coefficient  at  the  end  of  the  plate; 

Cf,L*=  0.455/(ln^(0.06p  VooL/ /)) 

=  0.455/(ln^(0.06(0.5528  kgW)(270.37  m/s)(l  m)  /  (2.091  x  10'^  N-  sW)) 

=  0.002705 

Eocal  Stanton  number  at  the  end  of  the  plate; 

StL*  =  (Cf,L*/2)  /  (l+12.7(Pr*2'^-l)(Cf,LV2)'^2 

=  (0.002705/2)  /  (l+12.7((0.6978)^^^-l)(0.002705/2)^^^)  =  0.001502 

Eocal  heat  transfer  coefficient  at  the  end  of  the  plate; 

/2l=  SfiVVooCp*  =  (0.001502)(0.5528  kg/m^)(270.37  m/s)(1009.44  J/kg-K) 

=  226.7  W/m^-K 

Average  heat  transfer  coefficient  over  the  length  of  the  plate; 

^  =  1.15  iL=  1.15(226.7  W/m^-K)  =  260.7  W/m^-K 

Average  heat  flux  dissipated  over  the  plate; 

q^=h  (Ew  -  Eaw)  =  (260.7  W/m^-K)(408.15  K  -  316.66  K)  =  23847.7  W/m^ 
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